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Chapter 1 


INTRODUCTION 


This thesis addresses the stability of coastal rubble mound stone armor 
layers exposed to water waves. The most common coastal rubble mounds are 
breakwaters, jetties, and revetments. Breakwaters usually have a large part of their 
length oriented perpendicular to the direction of wave travel and their primary purpose 
is to produce a quiescent area in their lee for ship anchorage or peed sheltering. Jetties 
are usually oriented parallel to the direction of wave travel and are used to produce an 
artificial channel or inlet. This eaaenel is usually an entrance to a harbor and is used for 
commercial or leisure water-craft navigation. Revetments typically armor coastal or 
riverside slopes that would otherwise erode when exposed to water waves. These waves 


could be a result of local or distant storms or they could be due to ship wakes. 


Breakwater failure can occur due to a number of different failure modes. 
The dominant failure modes are shown in Figure 1. Of these, seaside armor stability on 
a traditional rubble mound is critical to the integrity and functionality of breakwaters 
and is therefore the focus of the present study. Breakwater stone-armor layer stability is 


unique from other structural design in that it is highly variable, and this variability is 
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Figure 1.1. Breakwater failure modes 


complex and difficult to quantify for the wide range and many permutations of the many 
variables involved. The variability is dependent on the stochastic nature of both the 
armor boundary conditions and loading. Boundary conditions are the points of contact 
underneath and between the armor stones. The boundary conditions are uncertain 
because of the irregular stone shape and stone placement. The geotechnical 
characteristics of the foundation are also often uncertain. The loading of engineering 
interest is primarily due to incident storm waves that vary with storm intensity and 
storm location. The impact of the waves varies with water level and local bathymetry. 
The local water level varies with tide and storm surge. Waves can dislodge armor units 
by uplifting, rolling, or sliding individual units or by causing en masse movement of the 
entire armor layer. Armor units can be dislodged from the upslope layer without 


jeopardizing the integrity of the armor layer. But if enough armor units move and the 


underlayer is exposed, then the breakwater can erode quickly. This structural integrity 
threshold is uncertain because damage, defined by the eroded volume, may be focused 
at a point or distributed over a broad area. Because the toe forms a foundation for the 
armor layer, if the toe armor is mobile, then the entire armor layer can mobile. How the 


armor layer responds to armor movement at the toe is quite variable. 


The variability in both the loading and boundary conditions demands the 
consideration of the randomness of each stochastic parameter. Many clients, including 
the U.S. Army Corps of Engineers, now mandate that a risk analysis be done as part of 
major breakwater rehabilitation studies. This requirement is intended to provide a 
standardized engineering economic analysis technique for comparing competing 
alternatives. The engineering performance study that must form the basis of the risk 
analysis is often accomplished using a reliability analysis, where the reliability, or 
conversely, the probability of failure, is quantified for each alternative. In the reliability 
analysis, all alternatives that fall below a predefined level of reliability are either 
modified or discarded. For example, a breakwater may have several alternative armor 
layer designs including randomly-placed stone, pattern-placed stone, and randornly- 
placed concrete armor units. Each of these alternatives will have a different associated 
cost and a different probability of failure for the design single storm or sequence of 
storms. Determining the probability of failure eee instability of the armor layer 
requires knowledge of the rate of deterioration or the rate at which stones are displaced, 


which has not been quantified for breakwater armor layers. Reliability methods have 


only recently been adapted to breakwater design, and very little data exist to support 
prediction of rubble-mound deterioration. Despite the many hundreds, and perhaps 
thousands, of studies on breakwater armor stability, there have been few generalized 
studies of long-term deterioration of traditional rubble-mound armor layers. Moreover, 
there have been no generalized studies of deterioration due to variations in storm 
sequences using random waves. Finally, damage experiments to date have been 
primarily conducted with nonbreaking waves, which is atypical of many environments 


in which breakwaters are constructed. 


The purpose of this study was to first identify the primary mechanisms of 
stone instability and damage development when the breakwater armor layer is exposed 
to depth limited, breaking waves. Chapter 2 discusses initiation of stone movement. A 
relation for predicting the initiation of stone movement when exposed to vertical uplift 
of normally-incident plunging breaking waves is presented in Chapter 2. The equation 
is verified using data from a small-scale two-dimensional flume physical model study. 
This portion of the research effort qualitatively addresses incipient motion of stone 
armor. This section includes the traditional development of a stability prediction 


equation and some insight into the effect of wave steepness on instability. 


A primary goal of this study was to establish predictive relations for damage 
development on traditional breakwater sections for single storms and for storm 
sequences given depth-limited normally-incident waves. Chapter 3 introduces the 


subject with a discussion of historical stability and damage development physical model 


measurements and the relations derived from these past studies. In Chapter 4, a new 
flume study is discussed where breakwater profiles were measured as damage 
progressed on a breakwater cross section exposed to normally-incident depth-limited 
breaking waves. A new device for measuring breakwater profiles is discussed in some 
detail. New parameters are defined for prescribing the engineering characteristics of the 
eroded profile. In Chapter 5, measurements from this study are discussed in detail. In 
Chapter 6, profile and wave data are analyzed to produce spatial and temporal relations 
for predicting the mean and standard deviation of eroded area, eroded depth, eroded 
length, and remaining cover depth on a breakwater that is exposed to normally-incident 
depth-limited wave conditions. These equations should be valuable in support of 


reliability analyses as part of comprehensive risk analyses. 


Chapter 2 


INCIPIENT STABILITY OF BREAKWATER ARMOR UNITS 


2.1 Armor Stability Equations 

Extensive research on breakwater armor stability has produced many 
empirical stability models. PIANC (1976) provides a summary of early stability 
models based on regular wave experiments. PLANC (1992) provides a discussion of 
more recent irregular-wave-based models. The most widely known empirical stability 
model was developed by Hudson (1958, 1959), following the pioneering work of 
Iribarren (1938) and is typically seen in the following form 


y,H° 


2S De 
K,(S, - 1)° cot Gy) 


where 
W = weight of armor unit 
Y, = Specific weight of armor unit material 
H = design wave height at structure toe 


K, = _ tabulated empirical stability coefficient 
S, = specific gravity of armor unit material 


@ = seaside angle of armor slope relative to horizontal 


Hudson also expressed this equation in a slightly different form as 


H 
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where N, is the dimensionless stability number and D,,<) = (Ws)/y,)" is termed the 
nominal stone diameter. Hudson did not use the nominal stone diameter variable. It 
was used by both Thomson and Shuttler (1976) and Van der Meer (1988) who noted 
that D,;, is the length of a side of a cube with volume equivalent to that of the median in 
the stone weight distribution. A glance at equations 2.1 and 2.2 shows that stone 
movement due to wave forcing, characterized by the wave height, is resisted primarily 
by the stone weight. As shown by Hudson (1958), Ahrens and McCartney (1975), Van 
der Meer (1988) and others, the wave shape, structure porosity, armor shape, and 
structure cross-sectional shape all affect the wave force. Hudson also noted that the 
resistance to movement is affected by the friction between armor units, stone shape, 

_ upslope armor layer weight, and armor slope. Equation 2.1 is preferable from a physical 
perspective because it maintains the inertial form of the forces, which is appropriate for 
armor stability. But equation 2.2 combines the stability coefficient with the structure 
slope, which is desirable because the stability coefficient is a function of the structure 
slope (e.g. Shore Protection Manual (SPM) 1984). This second form also presents the 
stone stability equation in a form that is similar to sediment transport formulae. This 


will be discussed further later in this chapter. 


In equation 2.2, the stability number is a convenient scale for characterizing 
the incident wave height relative to armor stone size. It is typically in the range of 1 to 4 
for stable or marginally stable rubble mound breakwaters (Van der Meer 1988). The 
SPM (1984) provides Kp values for various armor layers at different slopes on 
breakwater trunks and heads exposed to breaking and nonbreaking waves. Zero-damage 
K, values are typically used corresponding to less than 2 percent by count, or five 
percent by volume, of the armor in a layer being displaced. Using this ngsdaunaee 
guidance from the SPM (1984), the stability number is in the range of 1 to 1.6 for 
angular stone armor layers at slopes of 1V:2H or steeper exposed to breaking waves. 
The SPM specifies K, values up to 2.2 for stone armor layers at slopes of 1V:3H 
exposed to nonbreaking waves. So for most stable coastal breakwaters, the stability 
number covers a narrow range from | to 2.5. For deformable structures where the 
armor stone is expected to be mobile, such as S-shaped breakwaters and berm 


breakwaters, Van der Meer (1988) suggests N, = 3 - 6. 


Hundreds of studies have been conducted to quantify the single empirical 
parameter K, in equation 2.1 for the wide variety of prototype conditions that might 
exist. Originally Hudson only explicitly included the effect of regular wave height, 
structure slope, and armor stone specific weight. Hudson found no clear effect of wave 
period on_armor stability for the nonbreaking regular wave conditions he and his 
colleagues tested. They simply determined K, corresponding to the lowest stability 


condition over a range of typical wave periods. The Hudson equation has been 


extended to include the effects of irregular breaking and nonbreaking waves and wave 
period (Ahrens 1975, Ahrens and McCartney 1975, Carver and Wright 1991, SPM 
1984), various armor layer types and armor gradation (e.g. SPM 1984), and number of 


waves (Medina and McDougal 1988). 


Other regular-wave stability formulations that have been utilized in recent 
years include Hedar (1960, 1986), Ahrens and McCartney (1975), and Losada and 
Gimémenz-Curto (1979). Ahrens (1975), Ahrens and McCartney (1975), Losada and 
Gimémenz-Curto (1979), and Pilarczyk and Den Boer (1983) all showed dependence of 
wave period on stability number for regular waves. Each of these authors showed that 
minimum stability occurred for surf similarity numbers or Iribarren numbers between 2 
and 4, corresponding to plunging to collapsing breakers. The surf similarity number 
was defined as & = tanO / (H/L)* where 0 = structure slope, H = regular wave height, 
and L = local or deep water wave length. Thompson and Shuttler (1975) conducted an 
extensive series of irregular-wave armor stability experiments. Their conclusions on 
damage progression were insightful and are discussed in the next chapter. Using 
Thompson and Shuttler’s data and data from his own experiments, Van der Meer (1988) 
developed a stability model which explicitly included wave period, structure 
permeability, storm duration and damage for a single design storm. The Hudson and 
Van der Meer equations will be discussed further in the following chapters. The above 
mentioned stability models predict minor damage reasonably well, although poor 


predictions are common. Pfeiffer (1991) compared the models of Hudson (1958), 


Hedar (1986), Losada and Gimémenz-Curto (1979), and Van der Meer (1988). Pfeiffer 
found that the Hudson equation and the Losada and Gimémenz-Curto equations 
performed the best when compared to field data, but none of the equations matched the 
prototype records all that well. Pfeiffer also noted that none of the stability models 


mentioned above are specifically suited for predicting extended damage. 


2.2 Armor Incipient Motion Studies 

The empirical stability models of Iribarren (1938), Hudson (1958), and 
others are based on a free body analysis of an armor unit undergoing forcing due to 
shallow-water waves. Early stability models assumed 1) the principal wave force was 
due to down- or up-rush on an unsheltered and unrestrained unit, 2) the wave force was 
drag dominant, and 3) the drag force would be critical if the maximum horizontal fluid 
velocity was used, which was considered to be proportional to the shallow-water 
incident wave celerity. But for an intact structure and prior to initiation of incipient 
motion, the armor units are typically partially hidden and restrained from up or down 
slope movement; so lift, inertia, and convection across the armor layer must be 
considered. Moreover, Sawaragi et al. (1982) showed that the maximum fluid velocity 
on a rubble mound was not necessarily proportional to the wave celerity. Sigurdsson 
(1962) made force measurements on armor composed of spheres set on extremely steep 
slopes with an impermeable underlayer and derived incipient equations of motion; but 
concluded by stating that the dominant mechanism of initiation of armor motion was 


still unknown and required further investigation. Although many authors have 
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measured wave kinematics on armor layers and the resulting forces, including Mizutani 
et al. (1992), Torum and Van Gent (1992), Torum (1994), and Cornett and Mansard 
(1994), there have been few observations of incipient movement of armor units 
discussed in the literature. As such, the relationships between incipient motion, wave 


kinematics, and forces on armor units are still unknown. 


Melby (1987) and McDougal et al. (1988) discussed a model for predicting 
the wave forces on dolos concrete armor units and the resulting incipient motion of a 
lone dolos in one of two orientations. Their model utilized linear wave theory and 
Morison forcing (Morison et al. 1950) with added mass coefficients for wave slamming 
from Kaplan and Silbert (1976) and Kaplan (1979). Kobayashi et al. (1990) presented a 
numerical model for predicting the displacement of armor on a traditional rubble 
mound. The shallow water wave model interacted with a permeable flow model and 
hydrodynamic drag, inertia, and lift forces were computed using a Morison formulation. 
_ The model was limited to forces parallel to the structure because only depth-averaged 
velocities were predicted by this one-dimensional model. Torum and Van Gent (1992) 
discussed a similar wave model and compared it to velocity measurements above a berm . 
breakwater. Torum (1994) discussed the measurements further. Although two- 
dimensional velocities were measured, vertical flow in the breaking wave was not 
modeled numerically. In addition, Torum noted that the inertial force was not well 
defined by the traditional inertia term of the Morison equation. Cornett and Mansard 


(1994) described an experiment where forces were measured on a panel of stones. This 


1] 


approach was unique and yielded insight into the average frictional force on sections of 
the armor layer. They found that the frictional force tending to dislodge armor units was 
greatest below the still water level and that the character of the force depended strongly 


on the type of wave breaking. 


This chapter discusses a series of physical model experiments to identify 
and develop predictive models for breakwater armor incipient motion and to relate this 
motion to existing empirical stability relationships. The experiments were conducted in 
wave flumes at the Waterways Experiment Station in Vicksburg, Mississippi. The first 
experiment consisted of measuring wave-induced fluid velocities on and within the 
armor layer and runup/down. In addition, free surface oscillations were measured while 
observing armor motion on stone and Core-Loc armor. The observations from this early 
study led to an incipient motion experiment using a fixed-sphere armor layer with 
several loose spheres placed at various depths within the armor layer. A dominant 


incipient armor motion mode and predictive stability equation were verified. 


2.3 Experimental Setup 

The initial experiments were conducted to determine the nature of armor 
incipient motion and surrounding flow. The instrumentation included a laser Doppler 
velocimeter (LDV), high-resolution video, and runup and vertical free-surface-piercing 
gages near and within the armor layer. The experiments discussed in this chapter for 


incipient motion were all carried out using regular monochromatic waves. Data analysis 


12 


was performed on short segments of between five and ten uniform waves to develop 


clear relationships between wave parameters and armor motion. 


The primary study was conducted in a 46-m long by 0.46-m wide by 1-m 
deep flume, with an offshore slope of 1V:30H (Figure 2.1). A conventional rubble 
mound cross section was constructed with various seaward slopes and armor types 
(Figure 2.2). Table 2.1 lists the different test plans. In this thesis, only the sphere and 
stone armor plans will be discussed (Test Plans 3 and 4) in order to Goetetn continuity 
with the stone layer damage discussion in subsequent chapters. For Plans 3 and 4, 
velocity measurements, sampled at 100 hz, were made throughout the water column 
from the toe to two armor dimensions above the still water level. The ranges of physical 


quantities and common dimensionless parameters for Plans 3 and 4 are listed in Table 
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0.46 cm 
ex. velocity 
measurement location 
Wave Board 


aa T | | _7 Runup Gages 


Array 2 © \ age 


| Flat 42.4m | 


Figure 2.1. Flume plan (top) and profile (bottom) views for incipient motion 
experiment 
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Table 2.2 
Ranges of Measured Physical Quantities and Common Dimensionless 
Parameters for Plans 3 and 4 


Parameter Range 


Wave height at toe, H, 2.4 to 18cm 
0.75 to 4.0 sec 
&-to 24 cm 


0.007 to 0.1 


Wave period, T 


Water depth at toe, h, 


Wave steepness, S,= H,/L, 
L, = deepwater wave length 


Surf similarity parameter, & = tan a/(H/L,)”* 0.1 to 1.4 
tan a = 1/30 = beach slope 


Surf similarity parameter, & = tan 6/(H/L,)* 1.5 to 21 


tan 6 = 1/2 = structure slope 
Relative depth, h/L, 0.009 to 0.28 


0.06 to 0.88 


Relative wave height, H/h, 


The LDV was a two-watt argon-ion two-component device assembled by 
the Dantec Corporation. Here, two-component describes the fact that the LDV 
measured velocities in two orthogonal directions. For this study, these directions were 
always vertically upward and horizontal, in the direction of wave propagation. The 
LDV worked in the back-scatter mode using a non-intrusive probe which contained both 
the emission and receiving optics. The benefits of this device included nonintrusive 
measurements, small measurement volume, clean drop-outs, high sampling rate, and no 
required calibration. Drop-outs were situations where there was no measurable 
backscattered signal (e.g. when the laser beams crossed at a point above water, after the 
wave crest had passed and the probe was above the level of the wave trough). During a 


drop-out, the LDV would produce a constant signal at the last measured value. Because 
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only the peak measured velocity values were used herein, drop-outs did not pose a 


particular problem. 


For this experiment, the LDV probe was pointed through the glass flume 
wall into the voids within the armor layer and measurements were made of the internal 
flow within the porous media. Many of the voids between armor units were more than a 
nominal armor diameter deep, so aaeagnnTENs could be made outside the flume wall 
boundary layer. One drawback to the LDV was the requirement of a full time operator 
with continuous attention to detail. Also, because of the small measurement volume, 
small changes in measurement location in the sheltered region behind a stone or within 
the armor layer often yielded large variations in measured average peak velocities. 
Therefore, the instrument required many measurements to map the flow field. So data 


analysis requirements were substantial for this experiment. 


The wave heights were determined using free surface measurements fie a 
vertical capacitance-type gage positioned at the location of the structure toe with no 
structure in place. Synthetic rubberized horse-hair mats were placed landward of the 
structure location to absorb the waves. The sampling rate for free surface measurements 
was 20 hz. The zero-downcrossing wave height was computed as the average height 


from a burst of approximately ten regular waves. 
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2.4 Incipient Motion Observations 

Several dominant incipient motion modes were identified during the stone 
stability experiment in Plan 3, Table 2.1. The following descriptions pertain to initial 
armor motion on an intact, as-built structure. Rolling was the only mode of motion for 
stones on the toe. Although both onshore and offshore motion was observed, the toe 
units always moved out of the layer in the onshore direction. Armor near the still water 
level was more likely to displace than armor in other areas. This appeared to be due to 
the fact that the armor was loosened in this area due to high velocities in the breaking 
wave jet. Once loosened, the motion would depend on the armor shape and its position. 
If the armor shape was flat, then the armor unit would flop back and forth until it rolled 
out of the armor layer, generally rolling upslope during uprush. If the armor shape was 
rounded, which was normally the case, the armor units would jump vertically under the 
steep wave face if the wave was severely plunging or collapsing. If the wave was 
surging, then loose units would only be displaced if they were exposed. There did not 
appear to be sufficient lift in downrush or uprush flows along the armor layer to displace 
the stones unless they were odd shaped (flat). The only displacement mechanism 
observed for rounded stones sufficiently hidden in the armor layer was uplift under the 


steep wave face. 


These observations indicated that a fluid velocity or acceleration component 
in the vertical direction is normally required to initiate armor motion for hidden armor 


units. Additionally, this early qualitative study indicated that, for a given wave height, 
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incipient motion was primarily influenced by wave shape, stone position, stone 


exposure, and stone looseness. 


2.5 Experimental Measurements 

Throughout the experiments, vertical and horizontal wave velocities were 
measured in the vicinity of the armor layer. Figure 2.3 ohms typical time series of the 
horizontal and vertical velocities on the structure measured for one run of plan 3 with 
the following characteristics: H, = 12 cm, T= 1s, d, = 24 cm, and d, = 8.8 cm, where H, 
is the average wave height measured at the toe, Tis the average wave period, d, is the 
toe water depth, and d, is the depth of the laser measurement. The laser measurement 
was made 1 cm outside the armor layer, rnedeived perpendicular to the outer armor layer 
profile line. The sign convention is such that the horizontal velocity is positive seaward 
while the vertical velocity is positive upward. Typically for these measurements, the 
horizontal velocity signal was considerably smoother than the vertical velocity signal, 
due primarily to the small amplitude of the vertical velocity relative to the horizontal 
velocity and the relatively large amount of turbulence near the armor layer. Figure 2.4 
shows a velocity vector time series over one wave period, measured 1 cm outside the 
armor layer profile line for he run of plan 3 with the following characteristics: H, = 8.4 
cm, 7=2 s, d,= 24cm, and d, = 17cm. Also shown is the wave profile at the point of 
maximum vertical velocity. The plot shows a large vertical velocity just below the steep 
wave front. Observed maximum stone movement for this wave profile position is also 


shown. Figure 2.5 shows an example of vertical velocities outside and inside the armor 
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layer. These measurements were also made during one run of plan 3 with the following 
characteristics: H, = 12 cm, T= 1s, d,= 24 cm, and d, = 8.8 cm for the external 
measurement and d, = 13 cm for the internal measurement. The measurement locations 
for these time series are shown in Figure 2.6, where the structure slope is 1V:2H. Here 
it is clear that the velocities within the armor layer are highly irregular due to 


turbulence. 


ah vertical velocity | 


-0.8 ny 


5 5) 6 6.5 7 US 
time (sec) 


Figure 2.3. Velocity time series for one wave period with H, = 12cm, T=1 s, d,=24 cm 
and d, = 8.8 cm measured 1 cm outside the armor layer 
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Figure 2.4. Velocity vector for one wave period for armor lifting measured 1 cm outside 
the armor layer with H, = 8.4 cm, T=2s, d,=24 cm, and d, = 17 cm 
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Figure 2.5. Vertical velocity time series outside (top) and inside (bottom) armor layer 
with H, = 12 cm, T= 1 s, d,= 24 cm, and d, = 8.8 cm for the outer measurement and d, 
= 13 cm for the inner measurement 
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Figure 2.6. Measurement locations for vertical velocities shown in Figure 2.5 


Figure 2.7 shows a typical plot of vertical velocity v,,/(gH,)”, measured 1 
cm outside the armor layer at various depths, versus the square root of wave steepness, 
where v,,, is the average of the highest one-third peak velocities for the burst of regular 
waves, g = gravitational acceleration, structure slope = 1V:2H, and L, = deep water 
wave length. For simplicity, v is used instead of v,,, in Figure 2.7 and hereafter. 
Relative laser depth, rd = d,/d,, is the ratio of the depth of the laser to the depth at the 
toe, measured from the still water level. As noted by Sawaragi et al. (1982), maximum 
non-dimensional velocities commonly occurred for collapsing to plunging breaking 
waves. Pilarczyk and Den Boer (1983) showed minimum stability for 1V:2H slope 


occurred for & = 3.3. In this case, this point occurs at (H/L,)” = 0.15, which is near the 
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maximum measured velocity for the waves with T= 2s. The peak vertical velocities for 


given wave period decreased with increase of rd in this figure. 


Vertical Velocity at Armor Surface 
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Figure 2.7. Maximum vertical velocity versus the square root of wave steepness for 
wave periods 7= | and 2 s and relative depths rd = 0.36, 0.5, and 0.7 
2.6 Incipient Motion Prediction 

The previous experimental results indicated that one of the dominant 
incipient motion modes was due to the vertical force occurring under the steep wave 
front. The balance of forces for vertical incipient armor motion with no external 
restraints yields the equality between the submerged armor weight and the vertical fluid 
force; W' = F,,. The vertical force at the steep wave front can be described by the 


Morison equation (Morison et al. 1950). 
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where p = fluid density, A = cross sectional area of armor unit in direction of flow, Cp = 
drag coefficient, v = peak vertical velocity just above the armor layer as plotted in 
Figure 2.7, V = armor unit volume, C,, = inertia coefficient, and dv/dt = total fluid 
acceleration. The drag force, the first term on the right side of Equation 2.3, can be 
expressed as a function of the armor nominal diameter, D,, by introducing an armor 


shape factor, K,, as follows 


A=K,D, (2.4) 


n 


The drag force in Equation 2.3 is then given by 


K,C 
ii, = pp Cupts Coe (2.5) 
where the nominal diameter was previously defined as 
13 
Di sV* = *) (2.6) 
i 


where W = armor weight and y,= armor specific weight. 


At the point of maximum vertical fluid velocity, the local vertical fluid 
acceleration, Ov/ot, and horizontal velocity, u, are negligible. As such, the total 


acceleration reduces to a convective term. 
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where x = horizontal coordinate and y = vertical coordinate. If we assume the 
convective acceleration across the armor layer to vary linearly vertically, then the 


acceleration can be expressed as 


dv eh OV K vie 
oh ‘ ae Tce (2.8) 


where K, is an empirical coefficient of order unity. The maximum inertial fluid force in 


Equation 2.3 can thus be reduced to 


F, = pD; C’,v? (2.9) 
with 
CHIC: (2.10) 


Substituting Equations 2.5 and 2.9 into the stability criterion W’ = F,, with W’ = pg(S,- 
1)D,7 yields a stability relation in form similar to Shields criterion (e.g. Raudkivi 1990) 


for the initiation of motion of sediment particles 


2 


< 


Vv 
so a (CoanGneye Dall 
D,g(S-1) (Crane) (2.11) 
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where S, = armor specific gravity, g = acceleration of gravity, and v. = critical vertical 
velocity at which armor just begins to lift. In terms of N, from Equation 2.2, Equation 


2.11 becomes 


Hi. gH, va -1 
eee To (CD + C;) (2.12) 


where H/, = critical wave height at toe. It is interesting to note that the stability number 
is primarily a function of the Froude number, v, / (gH,)*. This formula ties the 


traditional stability relations to local vertical velocity measurements. 


Based on results of detailed velocity measurements in the interior and just 
outside the armor layer, the vertical velocity gradient was found to be proportional to the 
ratio of the vertical velocity and the armor diameter, as assumed in Equation 2.8. The 
empirical convection coefficient is K. = 0.90 for this experiment. This is shown in 
Figure 2.8 for a group of experiments summarized in Table 2.3. All experiments listed 
in Table 2.3 were conducted with a seaward slope of 1V:2H, D, = 4.6 cm, and d, = 24 
cm. The velocity values are positive peaks from the aligned inner and outer vertical 
velocity time series. In Table 2.3, the velocity gradient Av/Ay = |(v, - v,/(y,- y,)|, where 
Vv, = outer peak velocity, v, = inner peak velocity, y, = inner velocity measurement 


elevation, y, = outer velocity measurement elevation. 
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Table 2.3 


Summary of Convection Measurement Experiment Results with cot 6=2, 
D,, = 4.6 cm, and d,= 24 cm. 


Measure- Outer 
ment Wave Velocity 
Location | Height | Period | Depth 
cm s 


near free 
surface 


one-third 
of depth 
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Av/Ay (s"') 


v,/D,, (5) 


Figure 2.8. Vertical variation of vertical velocity under a steep wave front across the 


armor layer as a function of wave steepness 


The drag and inertia coefficients can be more accurately defined if we 
assume spherical armor. Based on previous studies of forces on armor by Mizutani et 


al. (1992) and Torum (1994), reasonable estimates for drag and inertia coefficients are 


Cp = 0.8 and C,, = 0.4 yielding 


K,Cp 
Cp = 5 AUS (US, = zl) (2.13) 
Ci, = K.C,, = 0.36 (K, = 0.90) (2.14) 


m 


where K, = 1.21 corresponds to a sphere. The critical vertical velocity, v., for the 


incipient vertical armor movement reduces to 
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Ve 
———— = 1.2 (25) 
D, 8(5,-1) | 
where the critical vertical velocity, v., depends on the nominal diameter, D,, and the 
specific gravity, S,, only, for a loose armor unit. So, if the vertical velocity exceeds this 


critical velocity, motion of the sphere should occur. At the point of incipient motion, 


this critical condition can be expressed as 


y= v. = Dee (SIN) (2.16) 


Plan 4 in Table 2.1 was designed to test the above criterion. For Plan 4, the 
armor layer was constructed using silicon rubber spheres which were glued together and 
attached to an inflexible yet porous metal mat. The metal mat was placed directly on the 
underlayer and fixed to the flume walls. Several loose concrete spheres were placed in 
the armor layer along a line from above the still water level down to the toe. Each two 
loose spheres were separated by two glued spheres so that there was no interaction 
between loose spheres. The sphere layer of Plan 4 was constructed to have the 


minimum porosity of a sphere layer of 0.33. 


For Plan 4, the loose spheres would not move under any conditions unless 
they were slightly raised in the armor layer. This was accomplished by placing a 0.5- 


cm-thick spacer under each sphere. The primary effect of this was to raise the porosity 
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surrounding the loose sphere, providing a path for water motion under the sphere. The 
only motion observed for the raised loose spheres was vertical motion under the steep 
wave front, following a slightly elliptical path, and landing back in their hole after the 
wave front passed. For tests with vertical velocities corresponding to the critical value, 
the spheres were just lifting off. For the larger vertical velocities, the spheres were 
lifting entirely out of their initial holes, but settling back into their holes. This sphere 
motion under the breaking wave is shown in the sequence of photographs in Figure 2.9. 


The sphere in motion is just left of the black rectangle on the right side of the 


photograph. 


In the incipient motion experiment, spheres at a depth of one-third the toe 
depth were the most mobile while spheres at the still water level were somewhat less 
mobile. This movement corresponded to the variation of the vertical velocities in the 
water column as shown in Figure 2.6. Figure 2.10 shows the incipient motion criterion 
of Equation 2.16 versus wave steepness for Plan 4 using a few representative points 
from each motion category as summarized in Table 2.4. For this figure v. = 61.8 cm/s 
computed using Equation 2.16 with D, = 3.03 cm, S, = 2.083, and g = 980.6 cm/s. The 
dark horizontal line represents the theoretical incipient motion criteria while the velocity 
measurements are represented by the dark dots. Observed movement is noted for each 
data point. The vertical gaps between the lifting group of points and the stationary and 
rolling groups occurred because the vertical velocity increased dramatically under the 


steep breaking wave face. Therefore, it was difficult to get a continuous set of points 
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over the entire range of v/v. ratios. For the drag and inertia coefficients selected, the 


incipient motion criteria agrees reasonably well with the observed movement. 
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Figure 2.9. Photographs of sphere motion during a typical incipient motion experiment 
for Plan 4 
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L = sphere lifting 
FR = sphere rolling 
IS = sphere stationary 
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Figure 2.10. Incipient motion criterion (Equation 2.16 with v. = 61.8 cm/s) versus wave 
steepness 
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2.7 Conclusions from Incipient Motion Study 

Experiments on incipient motion of breakwater armor showed several 
modes of displacement. One dominant mode was due to vertical wave forces which are 
shown to occur at the point of maximum vertical velocity under the steep wave front. A 
simple relation was derived assuming a Morison-like wave force balanced by the armor 
unit submerged weight. The wave force model was composed of drag, due to the 
maximum vertical velocities, and inertia, due to the vertical convective accelerations. 
The maximum vertical convective acceleration is shown to be roughly linearly related to 
the square of the velocity, which puts the inertial force term into the same form as the 
drag term. The resulting incipient motion stability relation is similar in form to the 
Shields sediment motion criteria. Further, when expressed as a traditional stability 
number, incipient motion is shown to be a function of the Froude number, v/(gH)*. The 
incipient motion criterion shows promise in predicting the incipient motion of spheres 
for the conditions tested, but further experiments are required. This study provided 
some insight into how breaking waves can instigate armor motion and remove armor 
units from an intact armor layer. Further experiments were conducted using a high- 
precision force transducer to measure the forces on exposed spherical armor units. 


These measurements have not been analyzed completely and are not included herein. 
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Chapter 3 


HISTORICAL DAMAGE MEASUREMENT AND DESCRIPTION 


3.1 Damage Modeling Standards 

There is a substantial amount of literature concerning the measurement of 
damage on rubble-mound coastal structures. Hughes (1993) reviewed laboratory 
techniques for measuring damage. He noted three types of experiments for accumulated 
damage: (a) long-duration tests, (b) accumulated-storm-impacts tests, and (c) residual- 
stability tests. There is overlap among these three and few standards appear to exist for 
these types of laboratory studies. Jensen (1984) noted that model storm duration should 
generally be specified to provide the equivalent of 8 to 10 hr prototype. This is 
sufficient if equilibrium damage, where further waves cause no additional damage, 
occurs in this time; but if not, a subset of tests should be conducted to determine the 
ultimate damage level. Hughes stated that tests should be repeated at least two to four 
times to develop sufficient statistical certainty in the expected outcome with more 


extensive testing performed if the variance is large. 


There are two dominant methods for damage measurement: (a) visual: 
counting the number of individual armor units that have been dislodged and moved 


more than one nominal diameter from their original location and (b) profile 


58 


measurement: determining the eroded armor volume through profiling. A subset of the 
profiling method of characterizing damage, noted by Van der Meer (1988) for dynamic 
stability, is through description of the profile geometry. Torum et al. (1979) and Davies 
et al. (1994) described measurement of the minimum depth of cover, d. shown in Figure 


3.1, which is a reduction of the profile shape to a single parameter. 


ERODED AREA, A, 


d. 
DEPTH OF C 


Figure 3.1. Sketch of breakwater profile with definition of eroded area and depth of 
cover 
3.2 Damage Measurement Methods 
3.2.1 Eroded volume method 

The U.S. Army Engineer Waterways Experiment Station’s (WES’) 
historical method for characterizing damage utilized profiles to determine the 
percentage volume of stones eroded relative to the total volume of stones in the active 


armor layer. Hudson (1959) used this volume method. But the method apparently is 
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not described in detail in any public references. As such, it will be described herein. In 
this volume method, the active armor layer was defined as extending from the middle of 
the breakwater crest to one zero-damage wave height below the still water level. The 
damage was determined through profiles using a sounding rod with a circular foot of 
diameter equal to 0.56D,9, where D5) = (Wso/y,)'” is the nominal diameter of the 
median stone weight, Ws, and y, is the specific weight of armor stone. The sounding 
disc size was determined so that the before-testing armor layer thickness, determined 
using the measured profile, coincided with the theoretical value. The soundings were 
generally obtained on a horizontal grid spaced evenly at 1.5D,<). A number of profiles 
along the breakwater length were averaged to determine an average profile. The 
average damaged profile was subtracted from the undamaged average profile to get an 
average eroded area over the active region. The eroded cross-sectional area is defined in 
Figure 3.1. This eroded area was divided by the total area of armor in the undamaged 
average profile to get a percent damage D%. Hudson’s (1959) zero-damage criteria 
corresponded to D% < 1 percent. The zero-damage criteria given in the Shore 
Protection Manual (SPM 1984) corresponds to D% < 5 percent by the eroded volume 
method or 2 percent by count. The justification for the less restrictive zero-damage 


criteria is not clear but evolved over many years. 


The primary weakness of the WES eroded volume method is that, because 


the damage is only computed over the active region, the damage value will depend on 
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the structure geometry. This method also provides no indication of the severity of 


damage as characterized by the profile shape. 


Jackson (1968) gave damage due to regular waves in terms of the WES 
eroded volume definition of damage. Jackson’s damage values, given in Table 7-9 in 
the SPM, have been widely used to predict damage but are limited to regular waves with 
damage starting from an undamaged structure. Also, the damage is not given as a 


function of time, which is critical for determining the reliability of a structure. 


Broderick and Ahrens (1982) and Van der Meer (1988) defined a 


dimensionless damage index using profile data as 


Se peace = A. 
2/3 2 
Mg, Di s0 (3.1) 
Pa 
where 

A, = eroded volume per unit length or cross-sectional eroded area 
M,, = mass of median armor unit in mass distribution 
Pa = armor unit density 


where their variable nomenclature has been modified to avoid confusion. Broderick and 
Ahrens stated that A, was calculated from the profile data by determining the difference 
between before and after testing profiles. The difference in Broderick and Ahrens’ 


method from the traditional eroded volume method was that they nondimensionalized 
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by the square of the nominal stone diameter rather than the area in the undamaged armor 


profile. 


If the eroded area over the entire structure is used to compute S rather than 
just the active region, the ratio of damage by the eroded volume method to that of the 


damage index can be computed as per Cornett (1995) 


ee h a 

t —— 

So EO ee, (3.2) 
D% 100 D2, 


where 
t, = armor layer thickness 
w. = crest width 
h, = breakwater crest elevation above bottom 
h, = water depth at toe 


H = design wave height 


6 = seaside angle of armor slope relative to horizontal 


Equation 3.2 assumes the wave height is less than the depth at the toe. 
Cornett noted that the range of this ratio is 0.6 to 1.25 for typical rubble mounds. For 
S/D% = 0.8, he notes that the zero-damage criteria of D = 5% corresponds to0 < S <4. 


This is quite a broad criterion. The damage index method appears to give a better 
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representation of damage. Broderick and Ahrens and van der Meer noted a zero- 


damage level of S = 2. 


Thompson and Shuttler (1976) also used the eroded volume method. They 
used a profiler with a foot diameter of D,,<9/2 with soundings spaced at D,,<9, which is 
similar to that used in the WES experiments. The structure was surveyed after each 
1,000 waves. They computed the eroded volume V, using the trapezoidal rule and an 
average profile. A damage number N, was calculated, assuming a spherical armor 
shape, as the number of stones eroded in a 9D,,. wide region of the breakwater section 


as follows 


N Pate 
———— 
Tope (3.3) 
p,—Daso 
which is equivalent to 
54p, A, 
Ne PRM (3.4) 
™Pa D,s0 


over a 9D,. width of the structure, where 


p®,= armor bulk density 
V, = average eroded volume 
Pp, = actual armor density 


The difference between the method of Thompson and Shuttler and the 


damage index method is the 54p”,/mp, in Equation 3.4. For Thompson and Shuttler’s 
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tests, the density ratio was p,/p®, =1.81. Therefore, the method of Thompson and 
Shuttler yielded a damage index approximately nine times that of Broderick and Ahrens, 
or the width of their structure in nominal diameters. Thompson and Shuttler also 
determined the minimum armor layer thickness at failure. They defined failure as the 
point at which an area of exposed underlayer of diameter D,;. occurred. Note that the 
minimum armor layer thickness will not be zero at failure because it is expressed as a 


spatial average of several profiles. 


H. R. Wallingford, Ltd. (1990), showed that Equation 3.1 yielded very 
different results if a slightly different method was used to compute the average eroded 
area. The first method they used was that described for the WES eroded volume 
method, where an average profile was used to determine an average eroded area. The 
alternative method was to sum the eroded areas from all profiles in order to compute an 
average eroded area. The difference between the two methods ranged from 2 to 82 
- percent. In general, the difference decreased as the damage level increased. Note that 


most authors do not describe the method used to compute damage. 


All damage methods discussed above share a common weakness, namely 
they compute the average damage, which may be concentrated in one pocket or spread 
out over several areas. Also, none of the methods give any indication of the profile 
shape, or more specifically, the maximum depth of erosion, which is certainly an 


important parameter for a multilayer structure. 
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3.2.2 Stone count method 

Hedar (1960), Owen and Allsop (1983), Hughes (1993), and Davies et al. 
(1994) describe measurement of damage through stone counts. Besides visual counting 
during the test, photo overlays and digital image processing software can be used to 
determine damage by stone counts; but these methods are relatively complex and time 
consuming. Stone count methods are useful for determining very low damage values 
but become inaccurate if more than a few stones begin to move or if movement is due to 
sliding rather than dislodgement of individual stones. Stone count suffers from the 
same weaknesses as the eroded volume method, namely that the spatial concentration of 
damage is generally not specified and the mentinca depth of erosion is not computed. 


Stone count is also somewhat subjective. 


3.2.3 Recent damage CaSHTETTCnE methods 

Davies et al. (1994) provided a review of laboratory techniques for 
measuring breakwater profiles and methods for characterizing damage. They described 
the WES damage D% (SPM 1984) as a visual displaced stone count and made no 
reference to the WES eroded volume method. They introduced a technique to compare 
the stone count method of damage measurement with the damage index method of 
Broderick and Ahrens (1982). For damage measurement, they used a semiautomated 
profiler that measured profiles with a small spatial sampling interval by dragging a 
wheel over the structure face. During an experiment measuring damage on a riprap 


armor layer, they computed the damage index, apparently using an average profile, and 
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they noted that stone counts were more accurate when only a few stones moved, but the 
volume method improved accuracy of the damage measurement for advanced damage. 
They also noted that the depth of cover d, was useful in describing the degree of 
damage. They noted that, for an armor layer thickness of approximately 2D,,<), when 
d,=D, so, the underlayer was visible through a hole D,,., in size, and when d,=0, 


significant damage to the underlayer had occurred. 


3.3 Damage Measurement Experiments 

Historically, breakwater design has been accomplished using an empirical 
stability equation, such as the Hudson equation (Hudson 1958, 1959) as shown in 
Equation 2.1. As described earlier, for this equation, Kp is defined for a given level of 
performance, typically the no-damage condition represented by D% less than 2 percent 
by count or 5 percent by volume (SPM 1984). This technique assumed damage 
approached an equilibrium level of D%, where further regular waves at the design 
condition induced no further damage. This is based on regular wave experiments where 
damage reaches an equilibrium level or failure occurs relatively quickly. Hudson (1958) 
and Van der Meer (1988) expressed Equation 2.1 as a stability number as shown in 


Equation 2.2 or 


N, = (K,cot®)'? = 


(3.5) 


n50 


where A=S, -1. 
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The Hudson equation was defined for regular waves but has been extended 
to irregular waves using irregular wave physical model experiments. The SPM suggests 
Ho, the average height of the highest 1/10 waves, to replace the regular wave height in 
Equations 2.1 and 3.5, but the justification is oe Researchers attempted to relate 
regular and irregular wave effects on stability in the 1970s with no definite conclusions 
(e.g., Jensen (1984)). Recently, Vidal et al. (1991); Vidal et al. (1995); and Jensen et al. 
(1996) emphasized the need for characterization of the large waves in a random wave 
train. Vidal et al. (1995) noted that representation of the irregular time series by Hq, 
the average height of the highest 100 waves, provided a comparative level of damage to 
that produced by regular waves; but the equivalent Statistic of the form H,,, will depend 
on the number of waves. Medina and McDougal (1988) introduced an interesting, albeit 
not rigorous, method for interpreting Jackson’s (1968) regular wave damage results 
using a Rayleigh wave height distribution. Their method incorporated storm duration 
into the equation. In summary, the regular wave stability and damage formulations 
given above are conservative for design; but a universal analytical technique for 
extending these relations using irregular waves and resulting damage development has 


not been completed. 


Thompson and Shuttler (1976) performed both long-term deterioration and 
shorter single-storm damage tests using a riprap-armored embankment with an 
impermeable core. All of their tests were restricted to nonbreaking waves in front of the 


embankment, mostly deep water. The average run length was 5,000 waves, based on 
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mean zero-upcrossing wave period T,,, with intermediate surveys every 1,000 waves. 
They considered 5,000 waves to be a typical storm duration. Several of their significant 


conclusions concerning damage are paraphrased as follows: 


1. The rate of erosion of a given riprap is, as expected, 
strongly dependent upon “the significant wave height” 
H Wee 


2. The rate of erosion decreases with time and hence the 
damage history curves flatten out. At the lower wave 
heights, the curves can become nearly horizontal, giving 
an apparently stable riprap slope. 


3. The very long preliminary tests give no certainty of the 
riprap eroding to a totally stable equilibrium state, even 
with low damage rates. Thus, it is not safe to assume, as 
is often done in regular wave tests, that a slope will erode 
to stability. All that can be said is that the erosion rate 
may become small enough to be ignored in practice. 


4. The method of laying the riprap has a significant effect 
upon the damage history. 


5. The wave energy spectral shape as specified by the 
width parameter, €, does not affect the ultimate erosion 
volume. 


6. The failure criterion requiring a given area of exposed 
filter layer was easier to assess than that requiring the 
observation of the erosion of filter material and gave 
erosion volumes at failure which were independent of the 
filter grades used. 


7. Experimental time limited the maximum value of N,,, 
the average number of waves incident on the riprap, to 
5000, which is typical of a storm. In most cases this was 
too few waves to determine whether or not equilibrium 
damage was achieved or whether, at a given value of H,,,, 
the slope protection would eventually fail. 
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8. Even in the limiting case of waves with a small 

significant height incident for a very long time on 

relatively large riprap, there will be a few rare waves high 

enough to remove the smallest stones of the riprap pack 

and hence give damage. 

9. The movement of the stone is greatest on the flatter 

slopes although the net erosion is small. This movement 

results in self-healing by the smaller stones. 

These conclusions offer a somewhat different view of damage development 

than was accepted at that time based on regular wave experiments. The conclusion that 


an equilibrium level of damage may not occur provided motivation to include damage in 


a stability model. 


Using Thompson and Shuttler’s riprap stability data, Van der Meer (1988) 
stated that the damage rate should be linear up to 500 to 1,000 waves but “for large N,,, 
numbers a limit to the damage should be reached (equilibrium).” These criteria for the 
relation between damage and the number of waves were limited to tests where damage 
was larger than S = 3 after 5,000 waves and where the filter layer was not visible after 
5,000 waves. Van der Meer’s discussion of damage progression is limited to widely 
graded armor, which may not deteriorate in a manner similar to uniformly sized stone 


armor. This will be discussed further in the next chapter. 


3.4 Damage Progression Prediction 
Table 7-9 in the SPM (1984), included herein as Table 3.1, provides a 


deterioration model for armor stability based on regular wave data from model studies 


discussed in Jackson (1968). The table gives damage as a function of wave height 
relative to the zero-damage wave height for several armor types. The tabulated values 
can be used to formulate an equation for armor stone damage. As an aside, note that the 
dolos damage progression in the table is limited to stability and does not include 
breakage. But because dolos movement causes breakage, this damage progression may 


not be conservative. 


Table 3.1 
Jackson (1968) Breakwater Damage from 
Table 7.9 of SPM (1984 


Percent 
Damage 
D% 


anal eae 


40 to 50 


A simple empirical model for the best fit line through the data points in the 


table for rough stone damage is given by 


pa =| 


do 


| = 2Bex 2 us | - 18.08 (3.6) 
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The experiments supporting Table 3.1 were conducted with regular waves, so measured 
damage approached an equilibrium after a relatively short duration. Therefore damage 
is not a function of storm duration in Equation 3.6. Figure 3.2 shows data from Table 
3.1 and damage given by Equation 3.6. The relation utilizes H, the regular wave height, 
and H,,,, the no-damage regular wave height corresponding to 0 to 5 percent damage by 
the eroded volume method. This damage relation shows how damage varies with wave 
height but it assumes starting with an undamaged structure and damage approaching an 


equilibrium level. 
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Figure 3.2. Damage, characterized as eroded volume, as a function of monochromatic 
wave height for angular stone (SPM 1984) 


Van der Meer (1988) reanalyzed data from Thompson and Shuttler (1976) 
and conducted a number of additional experiments and found the damage index S to be 


related to the number of waves N,,,, the significant wave height H,,,, and the mean wave 
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period T,,. Van der Meer noted that S = 2 in Equation 3.1 provides a good estimate of 
the initiation of damage and that failure occurred when S = 8 for structure slopes of 
1V:1.5H and 1V:2.0H, where failure was defined as exposure of the underlayer through 
a hole D,,. in diameter. Van der Meer used test durations of 1,000 and 3,000 average 
wave periods. Virtually all of his tests were performed with nonbreaking waves. Van 
der Meer performed eight tests with depth limited waves and eight more where perhaps 
only the highest waves in the distribution were depth limited. Based on these data, he 
proposed breaking-wave-induced stability after N,, waves could be determined by the 


following equations. 


For plunging waves: 


H 0.2 
SeetS = 174\(62) P01 || eats (3.7) 
ADI, WN, 
and for surging waves: 
Eee Ss 0.2 Z 
—— = 14P°3) —_| cot8é, (3.8 
AD, so VN. | ) 
where 
Hyg, = wave height exceeded by 2 percent of waves in wave height 


probability distribution, where Van der Meer assumed H,,, = 1.4 H, 
to obtain Equations 3.7 and 3.8 from his formulations for 
nonbreaking waves at toe of structure 


A = (S,- 1) with S, = specific gravity of armor unit 
Ie = permeability coefficient 
6 = structure slope angle 
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a = tan0/(H/L,)” = surf similarity or Iribarren parameter 


H, = significant wave height which is equal to H,,, 
Ib, = piicn|(2 7) 
T = mean wave period 


Van der Meer’s breaking wave tests were performed with a traditional 
multi-layer rubble mound where the notional permeability was P = 0.5. This 
permeability is an empirical parameter without any regard to the flow throughout porous 
media. Using this value of P, Equations 3.7 and 3.8 suggest that S is approximately 
proportional to H, and N,,°°. Because the number of waves is defined as N,, = £/T,,, 
where f, is the total run time, these formulations indicate that damage increases with the 
square root of time. But because the Iribarren parameter is raised to a negative constant 
in Equation 3.7 but raised to a positive power of P in Equation 3.8, the effect of wave 


period is not clear. It is clear that damage progression is very sensitive to wave height. 


Equations 3.7 and 3.8 are of limited practical use for depth-limited breaking 
waves primarily because the supporting breaking wave experiment was extremely 
limited in scope. These equations are essentially the same as Van der Meer’s 
nonbreaking wave equations except the Rayleigh wave height distribution assumption 
of H, = H>,,/1.4 has been substituted. Van der Meer included the structure slope in the 
Iribarren or surf similarity parameter (Battjes 1974), but the beach slope is critical for 
depth-limited breaking waves. Van der Meer conducted experiments for only one beach 


slope and used a very narrow range of wave periods. Therefore, the effect of varying 
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beach slope or breaker type was not determined. For breaking waves, as the beach slope 
steepens, the wave breaking becomes more vigorous (Battjes 1974) and the slamming 
forces increase dramatically as the waves start plunging to collapsing (Bruun 1985). 
Van der Meer’s spilling breakers were not very severe with respect to stability relative 


to plunging or collapsing breakers. 


Equations 3.7 and 3.8 do not indicate decreasing stability as the depth to 
wave length ratio decreases and the wave breaking becomes more severe. Carver and 
Wright (1991) showed that stability decreases dramatically with decreasing relative 
depth. The minimum stability condition occurs where the wave breaking is 
characterized by plunging to collapsing breakers at the toe of the structure. Figure 3.3 
shows some of Carver and Wright’s data replotted with Hudson stability coefficient, Kp 
in Equation 2.1, as a function of relative depth. Here the water depth h, at the structure 
toe is normalized by the local wave length L, based on the peak spectral period 
~ computed using linear wave theory. It is clear that the Hudson stability coefficient 
decreases with decreasing relative depth. This is a reflection of the severity and location 
of the breaking wave and the resulting wave forces. As was shown in Chapter 2, the 
critical wave forcing for incipient motion was uplift occurring under the steep wave face 
for collapsing to plunging breakers. For depth-limited conditions, as the water depth 
decreases relative to the wave length, the wave face steepens, and vertical forces at the 


steep wave face are able to loosen and mobilize the stones. Therefore, for depth-limited 
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waves, the critical condition for stability occurs for relatively long waves that plunge on 


the structure. 


A minor difficulty with Equation 3.7 and 3.8 is a potential source of 
confusion because the stability number is given in terms of H,,,, but the Iribarren 
parameter is in terms of H,. Finally, all of the previous damage empirical equations 
listed above are limited to constant incident wave conditions and water level as well as 
the initial condition of an undamaged structure. These equations are intended to predict 
damage for a single design storm. Therefore, they cannot be used to predict damage 


development over the life of a structure as is required in a life-cycle cost analysis. 


3.5 Variability in Stability Results 

Carver and Wright (1991) conducted irregular wave stability experiments 
primarily intended to determine random variations of damage initiation due to varying 
armor placement and due to variations in wave time series realizations with constant 
spectral parameters. Their tests progressed to low damage levels, varying up to 7.7 
percent displacement, by count, of armor. The wave conditions were assumed to be 
constant for a given test. But they did vary the wave height and period between 
structure rebuilds if the damage had not progressed far enough. This testing strategy is 
not typical and their interpretation of the data assumed damage caused by one wave 
condition would be independent of the damage caused by the following wave condition. 
They made note of the uncertainty in stability of traditional stone-armored breakwaters. 


In Figure 3.3, for each value of relative depth, the stability coefficient is shown to 
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- scatter widely. They concluded that the uncertainty was due to differing construction 
techniques from test to test and different random number seeds being used to generate 


the spectra. 


Font (1968) noted that the armor placement technique is important for 
determining when initiation of damage occurs, but is less relevant for advanced damage, 
which is similar to conclusions stated above by Thompson and Shuttler. It is anticipated 
that the scatter in damage will be large for advanced damage primarily because of the 
random nature of the interaction between irregularly shaped stones. Melby and Mlakar 
(1997) showed how Carver and Wright’s data could be used to compute the reliability 
of a breakwater with respect to the zero-damage condition of armor stability and how 
Van der Meer’s relation could be used to compute reliability of a damaged mound. But 
little work has been done on computing the reliability of a damaged rubble mound 


through simulated damage using historical storm data. 
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Figure 3.3. Hudson stability coefficient versus relative depth for angular armor stone 
exposed to irregular waves (Carver and Wright 1991) 
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Chapter 4 


DETERIORATION EXPERIMENT 


4.1 Overview 
The deficiencies of previous damage experiments were noted in the 

preceding chapters. Because accurate quantification of damage is required for a risk 
analysis, a small-scale laboratory experiment was conducted to quantify breakwater 
deterioration. The initial experiment utilized a traditional trapezoidal multilayer rubble- 
mound armored with uniform-sized stone which was exposed to depth-limited irregular 
waves. The objectives of the experiment were as follows: 

a. Determine deterioration rates for the structure under long-term exposure to 

storm conditions. 
b. Develop and evaluate methods for damage characterization. 


c. Quantify development of damage with changing storm parameters, e.g., water 
level, wave height. 


d. Quantify profile development as related to damage. 
e. Quantify alongshore variability in damage. 


f. Evaluate and standardize methods for breakwater profiling. 
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4.2 Experimental Setup 
4.2.1 Flume and structure 

The experiment was accomplished using two identical side-by-side 
structures constructed in a flume measuring 61.1 m long, 1.52 m wide, and 2.0 m high 
(Figure 4.1). The structure cross sections are shown in Figures 4.2 and 4.3. Each 0.76- 
m-wide structure was constructed on a flat bottom. The offshore concrete beach slope 
was 1V:20H. This steepness was selected to realistically produce severely breaking 


waves in the vicinity of the structure toe. 


Wave Generator Nearshore Wave Test 
Offshore Gages Structures Wave 
“7 Wave Gages Absorber 


Test Section 
Viewing Window 


Figure 4.1 Flume profile for damage experiments 
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2 Core: passed #3 (0.67 cm), 2 
retained on #4 (0.47 cm) 2.9 cm 


Figure 4.2 Model structure cross section for damage experiments 
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Figure 4.3. Photograph of model structures before testing 


In armor stability experiments, breakwater crest heights are usually set high 
enough to prevent overtopping from influencing the stability, unless overtopping is the 
focus of the study or the breakwater is low crested. But the composite downslope 
weight of upslope armor units generally contributes to the stability of the armor layer 
and increases as the slope lengthens; although under some conditions, longer armor 
layers can have a tendency to buckle or slide. Just considering stone stability, for a 
given wave height, taller structures are typically more stable, with respect to main armor 
stability, than short structures. But, in the prototype, the cost of the breakwater is kept 
to a minimum by setting the breakwater crest height such that minor overtopping is 
allowed for a design wave condition. Thus, experimenting with non-overtopped 


structures that have high crest heights typically will not be conservative with respect to 
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damage development. For the experiment discussed herein, the structure crest height 
above the toe of 30.5 cm corresponded to an elevation of approximately one significant 
wave height above the still water level for the larger wave conditions tested. This is 
typical of prototype structure crest heights. This crest height was specified so that 
results from this experiment could be used for analysis of prototype breakwaters with 
higher crest heights because the damage measured with the shorter structure should be 
relatively greater than that of a taller structure. The effect of wave overtopping on 


stability will be discussed in more detail later in this chapter. 


4.2.2 Experimental test conditions 

The test conditions for the initial study are summarized in Table 4.1. Note 
that wave period variation and armor gradation variation experiments are not listed in 
this table, but were added, and are discussed in Chapter 7. The waves had relatively 
long periods in this initial experiment, as this produces the lowest stability (Figure 3.3). 
The depth-to-wavelength ratio h,/L, at the structure toe was approximately 0.07 for this 
initial experiment, where L, was computed using linear wave theory and the spectral 
peak period. This shallow-water breaking wave is also typical of design conditions on 


United States coastlines. 


4.2.3 Armor layer and underlayer 
The mass distributions for armor and underlayers for the initial test series 


are shown in Figure 4.4. The armor layer consisted of uniformly hand-sized stones with 
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Table 4.1 
Experimental Conditions for Initial Damage Experiment 


Waterdepthattogh, 14.9-15.8cm 
Se abi ee Semen 
(Stmicturetelopeveoton dhl UT eh ee 
[Structure crestheight above bottom, h, —|30.5em 
[Nominal armor stone diameter, Din = |3.64em 

Weegee OD 

aioe gremt eaat Wi plew yl 


Median armor stone mass, M,, 
Stone density, p, 


Ms, = 128 g, while the underlayer was sieve sized passing 1.59-cm and retained on 
1.27-cm sieves. The ratio of armor median mass to underlayer median mass was 

(M5)q / (Ms0),, = 25, which was large relative to that typically used in the prototype. The 
gradation of the armor layer for this initial series was very narrow with D,./D,, = 1.05, 

_ where D,, is the nominal diameter corresponding to the 85-percent exceedance 
probability in stone mass and D,, is the nominal stone diameter corresponding to the 15- 
percent exceedance probability. The gradation of the underlayer was wider with 

D;,/D,; = 1.44. The core material was sieve sized, passing No. 3 (0.67-cm) and retained 
on No. 4 (0.47-cm) sieves. While not impermeable, this core material was sized small 
enough to prevent strong flow within the core. Equations 3.7 and 3.8 indicate that lower 


porosity in the underlayer and core results in lower stability. Therefore, the underlayer 


21) 


and core were selected to be at the lower end of the range of material that is used in the 


prototype. 
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Figure 4.4. Armor and underlayer stone mass distributions for initial experiment 


The amount of material in the armor layer is determined using the packing 


density equation, given in the SPM as 


N, 
yi nk, (1 - P,/100) (y,/Weo)”* (4.1) 
where 
N, = number of armor units 
A = unit area of breakwater slope to be armored 
n = number of armor layer thicknesses 
k, = layer coefficient 
P, = armor layer porosity in percentage 


W.,) = median armor weight 
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It will be convenient to refer to damage in terms of the nominal stone 


diameter, so the packing density equation can be rearranged to be compatible as 


Np, 
We (4.2) 


ED 
2 
D n5O 
where @ is the packing density coefficient. The stone armor layer was placed in a 


traditional two-layer thickness at an average packing density coefficient of @ = 1.2, 


corresponding ton = 2,k,= 1.0, and P,= 40 percent as recommended in the SPM. 


The core was placed by shovel, troweled to grade, and washed in place to 
naturally pack tight. This process simulates the natural washing action of waves during 
the construction sequence. Similarly, the underlayer was dumped onto the sloping core 
and lightly troweled to grade. Each armor stone was placed individually by hand, 
simulating crane placement as closely as possible. The armor stone placement during 
the experiment followed WES guidelines for random hand placement in the laboratory 
as follows: 

a. Stones were lowered vertically into position. 

b. Stones were placed on the slope rather than dropped. 

c. Stones were placed so that they touched their neighbors. 
d. Stones could not be pushed into a hole. 


e. Particular orientation of stones or placement patterns were avoided. 
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4.2.4 Model-Prototype similitude 

Many authors have discussed similitude requirements for breakwater 
stability including Dai and Kamel (1969), Hudson et al. (1979), Sakakiyama and Kajima 
(1992), and Hughes (1993). As such, only a brief review will be provided herein. In 
order for damage measurements in the physical model to be scalable, similitude must be 
maintained, where the ratios of the dominant physical forces in the model are the same 
as they are in the prototype. Because there is no mathematical relationship that 
describes armor damage, dimensional and inspectional analysis and the resulting 
empirical relations are relied on to establish similitude requirements (Hudson et al. 


1979, Hughes 1993). 


The first requirement for similitude is that the model be geometrically and 
kinematically undistorted. Therefore, the size and shape of the breakwater, armor and 
underlayer stone, surface characteristics of the stone, and the height and length of the 
waves must all be undistorted. As described in Chapter 2, the wave forces on the armor 
units act to cause damage while the armor unit self-weight and inter-unit friction act to 
prevent armor movement. Also shown in Chapter 2, wave forces are typically 
decomposed into fluid drag, fluid inertia, and buoyancy. The two important force ratios 
that include these forces are the Froude number and the Reynolds number (e.g. Hudson 
et al. 1979). The Froude number is the square root of the ratio of the inertial and 
gravitational forces or Fr = u/(gl)* where u is the fluid velocity, g is the gravitational 


acceleration, and / is a characteristic dimension of an armor unit or wave height. Froude 
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similitude demands that Fr,,= Fr,, where m and p refer to model and prototype, 
respectively. The Reynolds number is the ratio of the inertial forces to viscous forces 
Re = ul/v, where v is the fluid kinematic viscosity. So, in addition to geometric 


similitude, complete similarity in damage requires Froude and Reynolds similitude or 


lip. Re, 
Tago Ne Ee 


In stability models, Froude and geometric similitude are typically maintained leading to 


the relationship between length and time scales of 


N, = N, = JN, (4.4) 


where N,, N,, and N, are the time, velocity, and length scales, respectively. Under 
Froude similitude, in order to achieve Reynolds similitude, the following relationship 


must be satisfied. 


i, = RE (4.5) 


where N, is the viscosity scale ratio. Therefore, it is not practically possible to attain 
both Reynolds and Froude similitude with the use of water in small-scale stability 
experiments. Water is always used in stability models, yielding some Reynolds scale 


effects. Thus the viscous forces are relatively stronger in the model than in the 


prototype. 
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In order to maintain reasonably close similarity between model and 
prototype, armor stability models are usually constructed at a scale large enough to 
avoid significant Reynolds scale effects by ensuring that the flow around the armor 
remains turbulent. Dai and Kamel (1969) suggested that Re = (gH)'”D,... /v > 3 x 104 
will prevent Reynolds scale effects. Their results were based on large- and small-scale 
flume tests using regular waves. Van der Meer (1988) stated that the lower limit of 


Reynolds numbers should be in the range 


H_ D 
Re = Y8=s" = 5 1x10*— 4x10" (4.6) 
Vv 
to prevent Reynolds scale effects. This guidance is presently accepted practice. In the 
experiment discussed herein, the Reynolds number range is 3.1 x 10° < Re < 4.0 x 10* 


which satisfies the requirements suggested by both Dai and Kamel and van der Meer. 


4.3 Wave Generation and Measurement in Initial Experiment 
4.3.1 Wave generation 

Waves were generated based on the Texel, Marsen, and Arsloe (TMA) 
spectrum (Bouws et al. 1985). The deterministic spectral amplitude and random phase 
method was used to synthesize the time series that was used to drive the piston-type 
wave board. Long (1986) describes the computer program used to synthesize the time 
series. The software creates a variable spectral bandwidth Aw, by the method of Goda 


(1970), which allows closer spectral lines near the peak frequency. For this method, a 
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sequence of random phases 6, are specified over the interval from 0 to 2m. Then the 


Fourier amplitudes are found from the relation 


c, = ¥25,(@,) Ao, (4.7) 


given a discrete target spectrum S, as a function of angular frequency w,. The final step 
is to determine the discrete time series by inverse Fourier transforming the amplitudes 
and phases. This method provides a very close match with the target peak frequency 


and the spectral shape. 


For these initial test series, two wave board drive signals were generated 
corresponding to two water depths. The design conditions for wave generation are 
listed in Table 4.2. Note that the gain is generator specific but provides a crude judge of 


the energy levels because gain is roughly linearly related to wave height. 


Table 4.2 
Wave Generation in Initial Damage 
Experiment 


Gain 
as percent of 
Depth maximum 


cm 
Se ee ae Pes 
[3 aes | ae eo 


Generator 
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The random number seed value for the random phase generation was 
different for the two depths. These two signals were scaled to achieve three wave 
heights for each water depth by adjusting the gain on the output amplifiers for the wave 
board command signal as indicated in Table 4.2. The resulting six time series of the 
command signals sent to the wave board are shown in Figure A.1 in Appendix A. Note 
that short segments at the beginning and end of each series have been cut off in the 
analysis. The spectra for these six signals are shown in Figures A.2 and A.3, for the two 
water depths. The wave board was commanded at a rate of 20 Hz, and ramps of 5 sec 
each were placed at the beginning and end of the command time series to gradually start 


and stop the wave generator. 


4.3.2 Wave measurement 

Waves were measured using capacitance-type wave gages. Two arrays of 
three gages, offshore and nearshore, were used as shown in Figure 4.1. Nearshore gage 
layout is chown in Figure 4.5. Considering the largest wave measured herein, this gage 
location approximately corresponds to the recommendation by Goda (1985) where the 
design breaking wave is determined 5H, seaward of the toe. This is the travel distance 
of large breaking waves. For the toe depth of h, = 11.9 cm, the depth at the most 
nearshore gage was 16.5 cm, and the depth at the offshore array was 112.7 cm. For h,= 
15.8 cm, the depth at the most nearshore gage was 20.4 cm and the depth at the offshore 


array was 116.5 cm. 


Landward Middle Seaward 
Gage Gage Gage 
(1) (2) (8) 


—fon a 0.305 m 


Figure 4.5. Flume profile near structure with nearshore wave gage layout 


The gages were calibrated using the following automated technique. Each 
gage was attached to a Jordan servo-controller and motor. The wave gage data 
acquisition, control, and signal conditioning system consisted of a VAX PDP-11 
microcomputer connected to 12-bit differential input analog-to-digital converters, 
digital-to-analog converters, and analog filters. A program on the PDP-11 sent signals 
- to the Jordan controllers sending each through 10 calibration stops twice, covering a 
range of motion in excess of the largest wave expected in the irregular time series. At 
each stop, the voltage from each gage was recorded. For each gage, a calibration curve 
was automatically tested for fit through the calibration points using linear, quadratic, or 
cubic spline regression. The calibration program automatically chose the simplest 
technique with less than 0.25 mm deviation from the regression curve, which 


corresponds to approximately 0.1-percent maximum deviation over the full calibration 
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range. This was the maximum error allowed and the error varied because the calibration 


range was different for each gage. 


Measured nearshore water surface oscillation time series are shown in 
Figures A.4 through A.6 for waves of three significant heights at h, = 11.9 cm. Figures 
A.7 through A.9 show wave gage time series for waves of three significant heights at h, 
= 15.8 cm. Incident and reflected waves were resolved using these time series from the 
three-gage arrays using the technique of Goda and Suzuki (1976). The technique was 
modified as per Kobayshi, Cox, and Wuranto (1990) to determine the incident and 
reflected wave spectra and time series at the shallowest gage. Time series parameters 
were determined using the zero-upcrossing method. At frequencies containing little 
wave energy, the reflection analysis technique can yield poor estimates due to low 
signal-to-noise ratios. These regions are indicated by low coherency in the cross 
correlation between gage pairs. For this study, the cutoff frequencies were established 
to maintain the coherence above 0.3 in the cross correlation between gage pairs. The 


analysis technique was not sensitive to coherence cutoffs higher than 0.3. 


The resolved incident and reflected spectra for time series from the 
nearshore and offshore arrays are shown in Figures A.10 through A.12 for h, = 11.9 cm 
and in Figures A.13 through A.15 for h,= 15.8 cm. Figures A.10 and A.13 show an 
increase in peak energy density due to wave shoaling between the offshore and 
nearshore gages and an increase in wave energy for low ( f < 0.35 Hz) and high (f > 1.0 


Hz) frequencies. These figures show that wave energy is being transferred due to 
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nonlinear wave interactions from the frequency band 0.35 -1.0 Hz to frequency bands 
f<0.35 Hz and f > 1.0 Hz. Figures A.11, A.12, A.14, and A.15, for the larger wave 
heights, show a similar transfer of energy but also show a decrease in peak wave energy 
due to wave breaking. It is clear from these figures that the wave energy decrease due to 
wave breaking increased dramatically as the generated wave energy increased. It is also 
interesting that the measured reflected wave energy was much higher nearshore than 
offshore. This trend is primarily because wave breaking and bottom friction reduce the 
incident wave energy as waves progress shoreward while the reflected waves undergo 
little energy change as they progress offshore. This is consistent with the measured 


cross-shore variations of wave reflection from beaches (Baquerizo et al. 1997). 


Tables 4.3 and 4.4 list nearshore incident-wave characteristics for the six 
wave trains, each of 15-min duration, measured with the structure in place. Table 4.3 
lists the spectral analysis parameters, and Table 4.4 lists the time series parameters, both 
- at the location of the shallowest gage that was 0.91 m seaward of the toe. The values 
were computed using the analysis method discussed above with the three nearshore 
gages. In Table 4.3, the mean period T7,, is defined by the relation 


J fSul fdf 
= eee (4.8) 


1Snl fat 
0 


Tm 


where S,(f) is the energy spectral density of the incident wave and fis the frequency. 


The spectral significant wave height is defined here as H,,, = 4m,'” with m,= zero 
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moment of the incident wave spectrum. The average reflection coefficient is defined as 
R= [(m,)/m,]'” with (m,), = zero moment of the reflected wave spectrum. The change 
of wave reflection with damage measurement was negligible in this experiment. This is 
consistent with the reflection measurements made by Smith et al. (1992) on a horizontal 
bottom in deeper water in front of a 1V:3H slope. The local wave length L,, and the 
deep water wave length L,,, = g7,,/27™ were both computed using linear wave theory and 
the mean period. The Iribarren parameter in Table 4.3 is defined as 


E, = tan o/(H,,,/L,,)'” and in Table 4.4 as €,,, = tan «/(H,/L,,,)"” where tan « = 1/20 is 


the nearshore beach slope and H, = H,,; in Table 4.4. 


Table 4.3 
Nearshore Wave Statistics for Incident Spectra and Related 
Parameters with Structure in Place 


Reflec- 
Peak tion 
Period Coeff. 

T, 


p 
sec 


[an eee 
p2] ue | wa | 170 | 248 | 060 | 0.067 | 0027 | 0.30 | 
pa | ie | 142 | 172 | 248 | 060 | 0.066 | 026 | 029 | 
pe | ise | isa | 168 | 250 | 061 | 0078 | 0036 | 026 | 
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Table 4.4 
Nearshore Wave Statistics for Incident Time Series and Related 
Parameters with Structure in Place 


[ase | 


For a Rayleigh distribution of wave heights, H,,,o/H, = 1.27 and H,,,/H, = 
1.40. Computing these ratios using the values from Table 4.4, H,,,)/H, = 1.19-1.26 and 
Hi,,/H, = 1.30-1.42, indicating the slight reduction of wave heights due to wave 


breaking. 


Traditionally, the incident waves are measured without the structure in place 
(Hughes 1993). For this study, the waves in the flume were measured both with and 
without the structure in place. For the condition with no structure in place, an array of 
three gages was placed with its centroid at the location of the structure toe for 
calibration purposes. The wave gage array was moved to the location shown in Figure 
4.5 when the structure was in place: Wave heights measured without the structure in 
place were smaller than those measured with the structure in place, and the difference in 


H,,,. tanged from 6 to 24 percent. The difference between the two sets of measurements 
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was due to the fact that the wave generation system did not allow absorption of reflected 
and re-reflected waves within the flume; so wave energy built up in the flume during 
testing. It should be noted that systems that absorb reflected energy are only partially 
effective: and so it is virtually impossible to prevent wave energy increase due to re- 


reflected waves in stability tests. 


The wave energy in this experiment reached an equilibrium state quickly 
during each run due to the relatively long wave period. A simple analysis was done to 
determine the rate of increase of wave energy. For this analysis, the wave time series 
were divided into eight segments, each of 4,096 points or 3.41 min, except the first 
segment which was 2,048 points. Both the first and second segments started at the 
400th data point or t = 20 s and the final segment ended at the 16,784th data point or t= 
14.0 min. The segments overlapped by 2,048 data points, and each segment represented 
approximately 120 mean wave periods, except the first segment, which was about 60 
wave periods long. For each segment, the incident H,,, was calculated using the 
technique described above. Figures 4.6 and 4.7 show the wave height variation for h, = 
11.9 cm and h, = 15.8 cm, respectively. In these figures, the ratio of the incident 
segment wave height H,,,(t,,) to the overall wave height H,,,(¢,) is shown as a function of 
normalized duration t,,/t,, where t,, is the duration to the center of the segment and f, = 
15 min is the total run duration. As can be seen from the figures, the wave height 


reached an equilibrium very early during each run. Based on these observations, the 
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wave heights measured with the structure in place were used, as they were more 


accurate than those measured without the structure in place. 


—=— Wave 4 —=— Wave5 —=— Wave6é 


Figure 4.7. Wave height variation throughout Waves 4, 5, and 6 at 15.8 cm toe depth 
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The wave reflection analysis as used assumed a flat sea bottom and that 
linear wave theory was valid. But the bottom of the flume sloped at 1V:20H and the 
waves were Clearly nonlinear. An analysis was done to determine the sensitivity of the 
wave analysis method to selection of gage location and therefore beach slope. Using the 
depth at the center of the gage array in the reflection analysis, rather than the depth at 
Gage 1, resulted in no difference for H,,,, 2-percent difference for H,, and 7-percent 
difference for H,,,. The difference in depths between Gage 1 and the array center was 
12 percent. H, and H,,, were more stable than H,., in this analysis. H, and H,,, for Gage 
1 closest to the toe were approximately the same for repeated tests. This is important, as 
the damage was shown to be approximately proportional to the fifth power of the wave 
height in the previous chapter, exaggerating any errors in wave height measurement. 
Therefore, the more stable values of H, and H,,, are relied on in the analysis that 
follows. Comparison of measured wave heights with breaker index curves, such as 
Figure 7-2 of the SPM, showed that the maximum H,,,,. would just break at the mean 
period in the depth of the shallowest wave gage, approximately 1 m from the structure 


toe. 


The Iribarren parameter or surf similarity parameter (Battjes 1974) is 


defined as 


H (4.9) 
L, 


where L, = g7T7/27 and tan 0 is either the structure slope or the beach slope, depending 
on where the wave breaks. The wave height H is usually either H, or H,,,, and the 
period is the mean period for armor stability (Van der Meer 1988). The Iribarren 
parameter can be used to characterize the type of wave breaking with — < 0.3 
corresponding to spilling waves, 0.3 < € < 2 to plunging waves, 2 < — < 3 to collapsing 
waves, and & 2 3 to surging waves. Van der Meer (1988) and others have used the 
structure slope in Equation 4.9. This is because they were primarily testing with 
nonbreaking waves on the beach. But for the experiment described herein, the wave 
heights were depth limited. So the beach slope was the critical parameter dictating the 
type of wave breaking near the structure and was therefore used to compute €. It is clear 
from Tables 4.3 and 4.4 that, for these tests, the most severe waves were primarily 
plunging. This corresponds to the design condition for most United States shoreline 
applications. The relative depths show that the wave condition selected is near the 


worst case for stability according to Figure 3.3. 


The stability numbers shown in Table 4.4 were in the range N, = 1.6-2.5. 
This is well within the range of NV, = 1-4 for conventional breakwaters suggested by Van 
der Meer (1988). On the other hand, stability coefficients in Table 4.4 based on Hj,,, are 
in the range 3.4-13.2 in comparison to Kp = 2.0 recommended in the SPM for breaking 
waves on rough angular stone armor placed on a trunk. This indicates that the 


breakwaters in this experiment will incur some damage. 


We 


4.4 Effect of Overtopping on Stability 
The effect of wave overtopping on armor stability was investigated by Van 


der Meer and Daemen (1994). They gave a reduction factor for stable armor size as 


1 


k, = for 0<R~<0.052 
1.25 -4.8R, # G2) 
with 
R S 2nH 
Re a aa eee S, = uy (4.11) 
P H,\ 2% y eT: 


where R- is the freeboard or structure crest height above the still water level. There is an 
odd mix of wave height statistics in Equation 4.11. Substituting the conditions in the 
experiment discussed herein yields R, of 0.063 - 0.076 for the toe depth of h, = 11.9 cm 
and 0.047 - 0.058 for h, = 15.8 cm. Therefore, the 11.9-cm depth is outside the range of 
applicability of Equations 4.10 and 4.11. Actually, wave overtopping was negligible in 
the shallow depth of h, = 11.9 cm. For h, = 15.8 cm, the range of stable armor size 
reduction coefficients is 0.98 < k, < 1.0. Therefore, wave overtopping had negligible 


effect on stability of the armor layer. 


4.5 Armor Profiler and Profile Measuring Technique 

A structure profiler similar in concept to that used by Davies et al. (1994) 
was constructed. The profiler design and construction are described in detail in 
Winkelman (1998), who participated in the planning and setup of the project. 


The profiler, shown in Figures 4.8 and 4.9, consisted of eight aluminum profiling arms, 
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laterally spaced at 5 cm, which pivoted about a point on a sliding carriage. A rotating 
spheroid with diameter approximately D,,,, at the end of each profiling arm followed the 
structure surface as the carriage passed over the structure, from just landward of the 


crest to just past the seaward toe. 


Lee Side 


Breakwater Crest 


PPPLEED 


76.2 cm 


Breakwater Toe 


Sea Side 


carriage 
©) 


‘otational 
potentiometer 


Figure 4.8. Plan (top) and profile (bottom) views of breakwater profiler 
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Figure 4.9. Photograph of breakwater profiler 
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Note that acoustic and laser profiling instruments were available as 
alternatives to the mechanical system used, but the signal returns off the irregular stone 
surfaces were poor resulting in unacceptably high data analysis requirements and low 
accuracy. Also, these backscatter-type sensors were generally an order of magnitude 
more costly. The spheroid diameter was chosen to coincide with that of Davies et al. 
(1994) so that results could be compared and was nearly twice the diameter of the 
sounding disc in the WES eroded volume method. In the present case, the sphere 
provides varying resolution of voids depending on the depth of penetration between 
stones and so corresponds to a range of sounding disc sizes. No attempt is made herein 
to relate the sphere size to the WES sounding disc size since there exist no damage 
progression data from the historical records that are directly comparable with data 
obtained for the present study. The profiler arm pivot point was at a fixed elevation just 
above the crest of the structure with the profiler arms extending back toward the 
structure. The angle of an arm at any point during the profiling process was determined 
through recorded voltages from rotational potentiometers attached to the pivot on the 
carriage. The location of the carriage was determined using a linear potentiometer fixed 


to the flume wall. 


It was required that the profiler arms be just heavy enough to continuously 
follow the structure without plowing, or moving stones. Another requirement was that 
the structure be profiled without lowering the water level between runs. Therefore, the 


hinge point, rod length, rod and ball weight, and counter weight were chosen to provide 


Tal 


an optimal balance over the length of the slope, where the profiling rods were out of 


water at the crest and partially submerged at the toe. 


The power input to, and analog output from, each of the potentiometers was 
controlled by a personal computer (PC)-based data acquisition system. Analog-to- 
digital data conversion and signal conditioning and recording were performed using an 
Optim MEGADAC 5000. The MEGADAC had 4 MB of internal memory and could 
sample up to 250 kHz from 20 channels. The heart of the MEGADAC consisted of 
AD682SH modules that included 12-bit differential input analog-to-digital converters 
and programmable Bessel filters. The MEGADAC also included power amplifiers to 
provide power to the potentiometers. The digital voltages were recorded to a hard disk 
within the PC at 20 Hz. The speed of the carriage over the structure varied, so this 
sampling rate produced a variable spatial sampling resolution. But, in general, the 
horizontal cross-shore sampling resolution was approximately one sample per 


millimeter. Optim, Inc. TCS software was utilized to contro] data acquisition. 


Prior to each day’s tests, the profiler was calibrated by first recording the 
rotational potentiometer voltages with the rods positioned horizontally on a jig that was 
set using high-resolution tiltometers. Then the rods were swung through a known angle 
and the voltage differences recorded. The resulting voltage-to-angle ratios provided a 
daily specific calibration coefficient for each potentiometer. Prior to each profile, the 
output voltage for each potentiometer was set to zero on the fixed bar just above the 


breakwater crest. 
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The elevations of each sphere bottom above the toe, z,, and horizontal 
distances from an arbitrary starting point, x,, were determined using the potentiometer 
measurements. The profiler sphere horizontal coordinate was computed as x, = x, + 
L(1-cos ), relative to x, = 0 with @ = 0, and the profiler sphere vertical elevation above 
bottom was computed as z, = z, - Lsin @ - D/2, where L = length of the profiler arm 
from pivot point to center of sphere, d = angle of the arm from horizontal, D = sphere 
diameter, x. = transverse distance the carriage moved, and Z, = elevation of the arm 
pivot point. Figure 4.8 shows a definition sketch for this transformation, where the 
initial carriage position x, = 0 corresponds to the initial sphere position x, = 0 with @ = 
0. The measured sphere elevations were averaged over a 0.5-cm interval in the cross- 
shore direction to remove small variations in elevations resulting from minor stone 
settlements that occurred during the tests. This is an essential step as this “noise” will 


cause a bias in the eroded area calculation. 


For each arm, the eroded area was calculated using a Simpson’s Rule 
numerical integration of the damaged profile below the undamaged profile. The eroded 
depth was determined as the maximum slope-normal distance between undamaged and 
damaged profiles. This slope-normal distance was calculated using a numerical 
technique beginning at each point on the eroded profile and stepping along a line normal 
to the slope of cot 0 = 2 until the undamaged profile was crossed. Then the point at 
which the crossing occurred was interpolated to obtain the slope-normal distance. The 


slope-normal minimum depth of cover at a given time t was calculated as the minimum 
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vertical difference between eroded profile and underlayer multiplied by the cosine of the 


structure slope angle or 


d(t) = min[z,(x,t) - z,(x,t=0)]cos 0 (4.12) 


where min indicates the minimum value with respect to x, z,(x,f) is the damaged profile 
elevation at location x and time ¢, and z,(x,t=0) is the underlayer elevation at x before 
testing began. The numerical interpolation technique used to compute the eroded depth 
was not required for the depth of cover as the underlayer was smooth relative to the 


profiler probe size in comparison to the undamaged armor profile. 
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Chapter 5 


DAMAGE MEASUREMENTS 


5.1 Overview 

In all, five test series were aaadneied with two structures in the flume within 
this initial experiment. The first test was not repeated. The other four series included 
two originals and two repeats; so three unique series were conducted and will be 
described in Chapters 5 and 6: Series A’, B’, and C’. As noted earlier, additional test 
series were added to investigate wave period and armor gradation and they are discussed 
in Chapter 7. Table 5.1 describes the test sequences for Series A‘, B’, and C’. Tables B1 


through B3 in Appendix B provide a more detailed summary of the Series A’, B’, and C’. 


Table 5.1 
Summary of Test Series for Initial Damage Experiment 
Order of Wave 
Test Test Conditions Water-Level Test raat 
Series Type (Tables 2 and 3) Order 
Deterioration to 1,2,3,4,5,6 Low - High 28.5 
failure 


i Storm ordering 1,2,3,5,6 Low - High ey 
4.5,6.2.3 | High-Low | 9.0 | 
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For Series A’, lasting 28.5 hr, the general experimental strategy was 
designed to provide an indication of the long-term deterioration of a structure. The 
structure was exposed to waves until failure, where the underlayer was visible through a 
hole in the armor layer at least D,5) in diameter. For this series, the experiment began at 
the lower water depth and progressed through each of the three wave heights listed in 
Tables 4.3 and 4.4, reaching an apparent equilibrium (visually observed) at each. Then 
the water level was raised, and each of the three wave conditions at the greater depth 
was run until an apparent equilibrium was achieved. The structures were exposed to a 
total of approkimaately, 60,000 waves for Series A’. This is equivalent to 100 to 250 hr of 
accumulated prototype storms for mean wave periods ranging from 7 to 15 sec. The 
number of waves was computed by dividing the run length by the mean period of the 


incident wave time series. 


For Series B' and C’, the testing strategy was intended to simulate damage 
from a sequence of individual storms. So for B’ and C’, each wave-water level condition 
was run for approximately 4,200 waves, representing roughly10-hr prototype storms. 
For B’ and C’, each sequence was run twice. And because there were two structures in 
the flume, this testing provided an original and three repeats to determine statistical 
variabilities. The total run length for B’ and C’ was 8.5 and 9 hr or approximately 
18,000 waves. Damage did not progress to failure for Series B’ and C’. In Series B’, the 
three wave heights at the low water level were run first (Waves 1, 2, 3), followed by the 


two highest wave heights at the high water level (Waves 5, 6). Series C’ was similar 
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except the high-water condition was run first followed by the low-water condition 
(Waves 4, 5, 6, 2, 3). Throughout the experiment, the waves were run in 15 min bursts, 
with the water completely settling between wave bursts. Profiles were measured after 


each two of these bursts, every 30 min of waves. 


5.2 Damage and Eroded Profile Parameters 

In Chapter 3, damage was characterized by the nondimensional eroded area 
S defined by Equation 3.1, and this parameter is used to describe damage herein. For a 
damage analysis, the shape of the eroded portion of the slope may be characterized by 
three parameters of primary engineering interest shown in Figure 5.1: cover depth d,, 
eroded depth d,, and eroded length /,. The cover depth, the minimum remaining depth 
of the cover layer along a profile, was shown by Torum et al. (1979) and Davies et al. 
(1994) to be a useful parameter to characterize the reserve capacity of the armor layer. 
For this experiment, the nondimensional cover depth was computed as the minimum 
depth of armor layer remaining for each measured profile after each 30 min of waves 


and normalized by the nominal stone diameter as 


C=— (5.1) 


where d. was computed using Equation 4.12 and C is the normalized minimum cover 


depth. 
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Underlayer 


Nondimensional 
Damage Parameters 


Figure 5.1. Sketch of breakwater profile with definition of damage parameters 


In this study, eroded depth and eroded length are introduced to describe the 
shape of the eroded area. They are parameters of primary engineering interest, 
providing a description of extent of damage normal to and along the slope, respectively. 
The eroded depth is defined as the maximum difference between undamaged and 


damaged profiles measured normal to the slope and is normalized as 


E=— (5.2) 


The eroded depth is useful because it gives an indication of progress toward failure, 
particularly if failure is defined by exposure of the underlayer. The sum of eroded depth 


and cover depth is approximately equal to the initial cover layer thickness t, 


2 @ 2 OC) Dx (5.3) 
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which is only approximate because the maximum d, and the minimum d, may not occur 


at the same location due to variability in the original layer thickness. 


The eroded length /, is defined herein as 


l=— (5.4) 


pe ee (5.5) 


Here the dimensional eroded length /, is not measured directly, but is a derived value. It 
is roughly equivalent to the surface length of the eroded region along a profile. As seen 
in Figure 5.1, the eroded length characterizes the along-profile length of the roughly 


triangular eroded area and describes the along-slope extent of damage. 


5.3 Probability Density Functions of Damage, Eroded Depth, and Cover Depth 
The number of damage profiles measured at a given time during damage 
progression is denoted by N, in the following, where N, = 16 for Series A’ and N, = 32 
for Series B’ and C’. The N, values for S, E, and C are analyzed statistically to reduce 
the number of measured values and separate the variability along the structure from the 
temporal variation associated with damage progression. First, the mean and standard 


deviation of the N, values of S, E, and C are calculated. The standard deviations 0, 0; 
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and 6, of S, E, and C indicate the variability of S, E, and C among the N, profiles. The 


means and standard deviations of S, E, and C vary with damage progression. 


The N, values of S, E, and C are normalized as 


S* = 50 COB = 5 €'es (5.6) 


where the means, indicated by overbars, and the standard deviations, o, of S°, E’, and C’ 
are zero and unity, respectively. As an example, the mean and the standard deviation 


for Series A’ after each 30 min of waves were computed as 


= 1 P 
Bink NS) \ vee 16 
res er Co 
N, 2 
0, = =e (S, - 5), N, = 16 (5.8) 


where S, is the individual measured value of S. The probability density functions P(S°), 
P(E’), and P(C’) are calculated using the N, values of S", E”, and C’, respectively, where 
the range between the maximum and minimum values is divided into eight bins of 
constant spacing. Admittedly, eight bins do not yield a fine resolution, but the number 


of data points for this estimation is limited to N, = 16 or 32. 


In summary, the N, damaged profiles measured every 0.5 hr yield one data 


set consisting of means and standard deviations of S, E, and C and the probability 
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density functions of S*, E*, and C*. The number of data sets is 57 for Series A’ lasting 
28.5 hr, 17 for Series B’ lasting 8.5 hr, and 18 for Series C’ lasting 9.0 hr. For the 
undamaged profile at the beginning of each test series, S = 0 and E = 0, which yield 

5, = 0, 6,, = 0, & = 0, and o,, = 0, where the subscript zero indicates initial values. For 
an undamaged structure, the probability density functions of S and E are the delta 
functions with a spike at S = 0 and E = 0, but S* and E* defined in Equation 5.6 cannot 
be computed because o,, = 0 and o,;,=0. On the other hand, the normalized cover 
depth C for the undamaged profile is a statistical variable related to the placement of the 
armor layer. The N, undamaged profiles are analyzed in the same way to obtain the 
mean C,, the standard deviation 6,,, and the probability density function P(C*) at the 


beginning of each test series. 


The probability density functions of S*, E*, and C* exhibit large scatters but 
do not indicate any specific variations with damage progression and among the three 
~ test series. All 92 data sets of P(S*) and P(E*) are plotted together in Figures 5.2 and 
5.3, respectively, whereas all 95 data sets of P(C*) are plotted in Figure 5.4. The 


standard normal distribution given by 


P(y) = exp(-0.5y2)//27 (5.9) 


with y= S*, E*, or C* is also shown in each figure for reference. Note that a fit to the 
distribution was not attempted because each data set consisted of too few points to 


determine the fit error. The data points in these figures scatter about the normal 
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distribution. The approximate ranges of the data points in these figures may be chosen 


as 


DY 2823 8 DI ZIFS2A 2 Dl BOS AS (5.10) 


which are in the range from -3 to 3. The corresponding ranges of S, E, and C can be 
found using Equations 5.6 and 5.10. This approximate procedure allows one to estimate 
the lower and upper limits of S, EZ, and C in terms of their means and standard 


deviations, which vary with damage progression. 


e SerilesA' co SeriesB' ~ Series C' —— Normal 


Figure 5.2. Probability density function for normalized damage S* for Séries A’, B’, and 
(ex 
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Figure 5.3. Probability density function for normalized eroded depth E* for Series A’, 
B’, and C’ 


° Series A’ > SeriesB' v Series C’ —— Normal 


Figure 5.4. Probability density function for normalized cover depth C* for Series A’, B’, 
and C’ 
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5.4 Temporal Damage Development 

The means and standard deviations of S, E, and C are used hereafter to 
examine damage progression. Figure 5.5 shows the structures prior to Series A’. 
Figures 5.6 through 5.8 show the measured profiles at the beginning, midway through, 
and at the end of Series A’. In Figure 5.6, one profile and the average of all eight 
profiles on one of the structures at the beginning of Series A’ are shown. In Figure 5.7, 
the average profile is shown along with a photo of the structures at the midway point of 
the series. This figure corresponds to an average damage level of $= 6.5. It may be 
seen from the photo that the top stones in the cover layer have been eroded sporadically 
over a wide region near the still water level. Figure 5.8 shows one profile and the 
average profile at the end of Series A’, where the underlayer is exposed through a hole of 
D,59 diameter. Figure 5.9 shows photos of the structures following Series A’. It is clear 


that most of the armor layer has been mobilized. 
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Figure 5.5. Photographs of undamaged structures prior to Series A’ 
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Figure 5.6. Profiles at beginning of Series A’, damage level of S=0 
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Profiles midway through Series A’, damage level of S'= 6.5 
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Gage 1 underlayer, damaged, and undamaged profiles 
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Figure 5.8. Profiles following completion of Series A’, damage level of $= 12.8 
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Figure 5.9. Photographs of structures following completion of Series A’, damage level 
of S= 12.8 


95 


Figure 5.10 shows the number of waves N,, versus the measured average 
damage S'for Series A‘ for N, = 16. Here, the cumulative number of waves N,, is 
obtained using the mean wave period T,, from Table 4.4 and duration of Waves 1-6 
listed in Table B1 in Appendix B. Also shown are the alongshore variability of damage 
as one standard deviation above and below S. For Series A’, waves were run at the 
initial depth until the structure appeared to stabilize, which occurred after approximately 
22,000 waves. At that point, the water depth was raised. As noted in Table 4.4, the 
stability number ranged from 1.6 to 2.2 for the shallower depth and from 1.7 to 2.5 for 
the deeper depth. In this figure, the significant wave heights H, and associated stability 
numbers N, are listed across the top of the figure, while the water depths at the toe are 


listed across the bottom of the figure. 


2D depth = 11.9 cm 
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Figure 5.10. Number of waves versus mean damage + one standard deviation for Series 
A ' 
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Figure 5.10 shows that the rate of damage increased with wave height, and 
damage became more scattered as the damage increased. Although the wave heights 
were reduced when the depth was increased approximately midway through the run, the 
damage continued to increase. This was due to the fact that, when the depth was 
increased, several stones higher on the structure (at the new still water level) were 
displaced down into the damage area, which was centered at the lower still water level. 
This can be seen in Figure 5.10 as small jumps in eroded area after the depth change. 


These small increments of damage approximately correspond to one stone displaced. 


Although the structure visually appeared to stabilize at various times during 
Series A’, Figure 5.10 indicates that the structure never really stabilized or reached 
equilibrium, where further waves produced no additional damage. Thompson and 
Shuttler (1976) noted this for their riprap tests. It is also clear that an extraordinarily 
large number of waves was required to induce failure. For Series A’, approximately 
~ 60,000 waves at significant wave heights appreciably greater than the no-damage wave 
height were required to fail the structure. Again, failure was defined as exposure of the 
underlayer through a hole at least D,,, in diameter. A higher mean damage value at 
failure was recorded for Series A' than was noted by Van der Meer (1988). Van der 
Meer quoted a failure damage value of = 8, whereas = 13 at failure for Series A’. It 
was noted during the series that the waves had a difficult time moving the second layer 
of stones. This led to high damage values at failure, with most of the first layer of 


stones moved before the second layer began moving. Most other recent damage testing 
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has been done with riprap, which consists of stone having a wide mass distribution 
(Thompson and Shuttler 1976; Davies et al. 1994). Riprap may not exhibit this 
secondary armoring because the small stones in the mass distribution are continuously 


moving. Deterioration of a riprap layer is examined in Chapter 7. 


The sawtooth nature of the damage curves was caused by partial self- 
healing. If a stone moved out of its location and started downslope, the next profile 
would show a hole and relatively increased damage. But during the following wave 
burst, the hole would fill in through adjustment of surrounding stone so that the next 
profile would show a partially healed structure. This healing process can be clearly seen 
in Figure 5.10 just after the water level was raised. Over time, the structure became 
more and more loosened, to the point that it could no longer heal itself because the 


surrounding stones were also displaced. 


Series B' and C' were similar to Series A’ in that the wave and water level 
sequences were similar; but the durations and ordering of the storms were different for 
B' and C’. The results of Series B’ and C’ are tabulated in Appendix B. Mean damage 
versus number of waves for the storm sequences of Series B’ and C’ are shown in 
Figures 5.11 and 5.12. Also shown are the variability as a standard deviation above and 
below the mean. For these series, the 32 profiles were used to compute each data point, 
corresponding to the four repeats with eight profiles each. The characteristics of armor 
deterioration for Series B’ were similar to Series A’ because the order of the storms was 


the same for both series. But for Series C’, Waves 4-6 at the high-water level increased 
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Figure 5.11. Number of waves versus mean damage + one standard deviation for Series 
B' 
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Figure 5.12. Number of waves versus mean damage + one standard deviation for Series 
(ey 
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damage initially. Waves 2 and 3 at the lower water followed but caused little additional 
damage to the segment of the armor layer damaged by Waves 4-6. This suggests that 
the damaged armor layer can withstand subsequent wave action at a lower intensity if 


the wave action is confined to the damaged area. 


As one would expect, the damage rates differed for the different storm 
ordering. But the ultimate damage at the end of the storm sequence was surprisingly 
consistent. Both the mean damage and the standard deviation for the four repeats were 
consistent. At the end of Series B’ and C’ after approximately 18,000 waves, S'= 8 and 
0, = 2 for both series. Although the testing was by no means comprehensive, it suggests 
that damage and variability are similar after being exposed to similar cumulative wave 
action with different sequences. This is useful in determining the reliability of a 
structure because the cumulative damage may be assumed to be the sum of damages 


caused by individual storms. 


5.5 Characteristics of Profile Erosion 

For Series A’, the mean minimum cover depth was computed as the average 
of the minima from the 16 profiles, while the mean eroded depth was computed as the 
mean of the maxima from the 16 profiles. Figure 5.13 shows the mean eroded depth 
versus number of waves, and Figure 5.14 shows the mean cover depth versus number of 
waves, both for Series A’. These figures clearly show the periodic healing process as a 
jagged sawtooth shape to the curves. They also show that the armor layer thickness 


decayed to an equilibrium level, with an asymptote at C=0. Figure 5.15 shows that the 
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eroded length jumps to a level of L = 7 initially and then increases with the number of 
waves in a manner similar to damage, as one would expect. The initially large value of 
L computed from a small value of £ illustrates that there was quite a lot of stone 


settlement but not displacement as damage was initiated. 
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Figure 5.13. Number of waves versus mean maximum eroded depth + one standard 
deviation for Series A’ 
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Figure 5.14. Number of waves versus mean minimum cover depth + one standard 
deviation for Series A’ 
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Figure 5.15. Number of waves versus mean maximum eroded length for Series A’ 
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Figure 5.16 shows mean eroded depth and mean cover depth versus mean 
damage for Series A’. As one would expect, cover depth is inversely related to damage, 
while eroded depth is directly related to damage. Eroded depth and cover depth are 
clearly inversely related. Similarly, Figure 5.17 shows that mean eroded length and 


mean damage are directly related. 
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Figure 5.16. Mean damage versus mean maximum eroded depth and mean minimum 
cover depth for Series A’ 
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- Figure 5.17. Mean damage versus mean maximum eroded length for Series A’ 


Figures 5.18 through 5.23 show that the profile development characteristics 
for Series B’ and C’ are similar to those of Series A’ discussed above. Figures 5.24 
through 5.27 show means of the three profile parameters versus mean damage for Series 


B' and C’. It is clear that profile development is related to mean damage. This will be 


investigated in the following chapter. 
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Figure 5.18. Number of waves versus mean maximum eroded depth + one standard 
deviation for Series B’ 
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Figure 5.19. Number of waves versus mean maximum eroded depth + one standard 
deviation for Series C’ 
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Figure 5.20. Number of waves versus mean minimum cover depth + one standard 
deviation for Series B’ 
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Figure 5.21. Number of waves versus mean minimum cover depth + one standard 
deviation for Series C’ 
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Figure 5.22. Number of waves versus mean maximum eroded length for Series B’ 
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Figure 5.23. Number of waves versus mean maximum eroded length for Series C’ 
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Figure 5.24. Mean damage versus mean maximum eroded depth and mean minimum 
cover depth for Series B’ 
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Figure 5.25. Mean damage versus mean maximum eroded depth and mean minimum 
cover depth for Series C’ 
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Figure 5.27. Mean damage versus mean maximum eroded length for Series C’ 
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Chapter 6 


DAMAGE AND ERODED PROFILE PREDICTION 


6.1 Damage Variability 

Figures 5.10 - 5.12 indicated that the standard deviation of damage S 
varying along the structure increases with mean damage S. Figure 6.1 shows the 
relationship between o, and S for all data in the three series plotted in Figures 5.10 - 


5.12. The scattered data points for the three series may be represented by the relation 
o, =05S°” (6.1) 


which is the solid curve shown in Figure 6.1. This relation implies that damage 
variability along the structure increases with the mean damage. On the other hand, 
relative variability defined as the coefficient of variation V, = o,/.S decreases with 
increasing §. This implies that higher damage levels may be estimated with smaller 
relative errors. The correlation coefficients for the three series for Equation 6.1 were all 
r=0.99. An example of the use of Equation 6.1 is illustrated for S ~13 at failure of 
Series A‘. Substituting S= 13 into Equation 6.1 yields 0, = 2.65. The approximate 
range of -2.7 < S* < 3 in Equation 5.10 corresponds to 6 < § < 21, indicating damage 


variability of a factor of 3.5 along the structure. 
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Figure 6.1. Prediction of damage variability, characterized by the standard deviation, as 
a function of mean damage for Series A’, B’, and C’ 
6.2 Eroded Profile Prediction 

It is useful to express the normalized eroded depth E and the variability of 
- eroded depth, o;, as a function of mean damage also. The measured values of and F+ 
6, were shown as a function of cumulative number of waves for Series A’ in Figure 
5.13. The trends for Series B’ and C' were shown to be similar in Figures 5.18 and 5.19. 
Figures 6.2 and 6.3 show and a; as a function of S for all data points of the three 


series. These figures indicate the following empirical relationships 


E =0.46./5 C2) 


= 0.26 — 0.00007(S —7.8)4 (6.3) 


iy 
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Figure 6.2. Prediction of mean maximum eroded depth as a function of mean damage 
for Series A’, B’, and C’ 


0.4 
eee 
Savill bal oe 
203 
@® 
o aes | 
ao Vv be 
> 
() 
ot 
zg 0.1 4 
7) 

(0) 


0 2 4 6 Bo 12 14 
Mean Damage, S 


Figure 6.3. Prediction of standard deviation of maximum eroded depth as a function of 
mean damage for Series A’, B’, and C’ 
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The correlation coefficients for Equation 6.2 were r = 0.99 for all three series. Equation 
6.2 indicates E/.S°° is approximately constant. This implies that the average shape of 
the eroded area remains approximately constant during damage progression. Figure 6.3 
and Equation 6.3 suggest that the variability of eroded depth increases rapidly with 
damage and becomes approximately constant at = 4. The variability in eroded depth 


decreases somewhat as the damage approaches the failure level of the armor layer. 


The eroded length L along the slope depicted in Figures 5.15, 5.17, 5.22, 
5.23, 5.26, and 5.27 is analyzed in the following. The mean eroded length J is simply 
defined as { = 2A, /d,, corresponding to a triangular shape. The normalized eroded 
length L = £/D,,., is then given by L = 2.5/E, as shown in Equation 5.5. From the 
measured values of and E, the corresponding values of L are calculated to examine 
the variation of L during damage progression. Figure 6.4 shows all calculated values of 
L for the three series against the corresponding values of §. Figure 6.4 also shows the 


following relationship derived from Equations 6.2 and 5.5 


L=44V5 - (Se) 


The correlation coefficients for Equation 6.4 were r = 0.98 for Series A’ and r = 0.99 for 
Series B’ and C’. As an example of the use of Equation 6.4, for S' = 13 at failure for 
Series A’, Equations 6.2 and 6.4 yield E = 1.7 and L = 16, implying the average 
damaged profile at failure extends 16 stone diameters along the slope. For the model 


structure with D,,. = 3.64 cm shown in Figures 4.2 and 4.3, this damage extent 


WN) 


corresponds to most of the seaward slope and crest of the structure as shown in Figures 


5.8 and 5.9. 
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Figure 6.4. Prediction of mean eroded length as a function of mean damage for Series 
A’, B’, and C' 

Finally, the normalized minimum cover depth C is analyzed to examine the 
decrease of protection with damage progression. The measured values of Cand C+0, 
were plotted as a function of cumulative number of waves in Figures 5.14, 5.20, and 
5.21 for Series A’, B’, and C’, respectively. These relations are similar to the 
corresponding results for damage and eroded depth except that C decreases with N,, and 
Gc is positive at N,, =0. The initial values of C at N,, = 0 were C, =1.65, 1.52, and 1.38 
for Series A’, B’, and C’, respectively, and the initial value of 0, at N,, = 0 was 0,, = 


0.14, 0.17, and 0.22 for Series A’, B’, and C’, respectively. To account for the 
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differences of C, and 6,,, (Cy - C) and (G¢ - Gg,) are plotted against .§ for all the data 
points of the three series in Figures 6.5 and 6.6, respectively. The empirical 


relationships shown in these figures are given by 


C=C, 05 (Ge) 


(6.6) 


Mean Damage, S 


e SerilesA’ © SeriesB’ vw Series C’' —— Eq. 6.5 


Figure 6.5. Prediction of mean minimum cover depth as a function of mean damage for 
Series A’, B', and C’ 
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Figure 6.6. Prediction of standard deviation of minimum cover depth as a function of 
mean damage for Series A’, B’, and C’ 

For Equation 6.5, the correlation coefficients were r = 0.99 for all three 
series. Figure 6.5 and Equation 6.5 indicate that the mean cover depth decreases 
approximately linearly with mean damage increase. The scattered values of (G¢ - Gg,) 
plotted in Figure 6.6 are of the same order of magnitude as the initial value 6, listed 
above. Figure 6.6 and Equation 6.6 indicate that the standard deviation 0, representing 
the variability of cover depth along the structure increases somewhat with damage 
progression and then decreases as the cover depth approaches zero. It is noted that 


Equation 6.6 overestimates 6- somewhat for Series C’ as can be seen in Figure 6.6. 


A failure criterion based on the normalized minimum cover depth C is 


explored in the following. Using the lower limit of C* = (C - C) /o, in Equation 5.10, 
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the lower limit of C may be estimated as (C -2.70,). The failure may be assumed to 
initiate when this lower limit becomes zero. This criterion with Equations 6.5 and 6.6 


yields the following equation for the mean damage at failure initiation 


0.18 + 2.7[0.098 - 0.002(5-7)"] = C, - 2.70¢, (6.7) 
This criterion accounts for the initial minimum cover depth and its variability along the 
structure where (C, -2.76,,) = 1.27, 1.06, and 0.79 for Series A’, B’, and C’. Solving 
Equation 6.7 for mean damage for Series A’ yields S= 10.9 at failure initiation, 
compared with S$ ~ 13 measured. The failure initiation criterion of Equation 6.7 based 
on the lower limit of C* of all data points plotted in Figure 5.4 yields the lower limit of 
the mean damage at failure. The upper limit of mean damage at failure may simply be 


estimated by C = 0 in Equation 6.5. This upper limit is S= 16.5 for Series A’. 


6.3 Temporal Damage Development 

If the mean damage can be predicted, then Equations 6.1 - 6.6 yield o,, E, 
o,, L, C, and o,. These relationships have been obtained by analyzing the measured 
profiles of the armor layer and underlayer without regard to the incident waves and 
water level. This statistical analysis indicates that the damaged profile statistics can be 
represented by the mean damage S. The next step is to predict the temporal variation of 
Son the rubble-mound breakwater exposed to depth-limited breaking waves in 


sequences of storms with varying wave conditions and water levels. 


7) 


The empirical formulas of Van der Meer (1988) indicate that breakwater 
damage due to depth-limited waves is a function of N,, T,, and time t. In addition, it is 
likely that damage is a function of structure slope, structure permeability, armor 
gradation, armor porosity, armor stone shape, and method of armor placement. None of 
these latter parameters except armor gradation were varied in this study. Gradation is 


discussed in Chapter 7. The remainder of the parameters will be investigated in future 


studies. 


In the following, the depth-limited breaking-wave stability data of Van der 
Meer (1988) will be used to develop relations for damage progression. His depth- 
limited wave and measured damage data are shown in Table 6.1. The particular data set 
shown in Table 6.1 corresponded to a structure slope of 1V:2H, a beach slope of 
1V:30H, and the stone armor characterized by D,, / D,;= 1.25, D,5)= 3.60 cm, and 
A=1.615. The table shows that, for toe depth h, = 0.4 m, the waves were only 
marginally depth limited because H,/h, was in the range 0.26-0.31. Therefore, Van der 
Meer’s data provide only eight tests where the waves were clearly depth limited (h, = 
0.2 m). The temporal development of damage cannot be deduced from these data 


because the darnage was only given after 1,000 waves and 3,000 waves. 


From Van der Meer’s data, it can be seen that the surf similarity parameter 
with respect to the beach slope for the eight breaking wave tests 9-16 ranged from 0.20 
to 0.25. This suggests that all were spilling breakers, which is not the worst case for 


armor stability. As explained in the previous chapters, plunging breakers expose the 
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armor units to the highest vertical forces at the steep front of a plunging wave, which is 


the worst case for armor stability. 


Table 6.1 
Depth-Limited Wave Damage Data from Van der Meer (1988 
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Consider Equations 3.7 and 3.8, which are Van der Meer’s equations for 
stability that include damage as a parameter. The stability number was a minimum at 


the transition between plunging and surging waves occurring at 
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Enc = (6.2 P°?!Vtan8]? (6.8) 


where Q = (P + 0.5)' and P = permeability coefficient varying in the range 0.1 - 0.6. 
Keep in mind that 0 in Equation 6.8 is the structure slope measured from horizontal. 
Later in this chapter, this requirement will be abandoned. Using Equations 3.7 and 6.8, 


the minimum stability number can be expressed as 


= \02 
= S 
N, = C.. 6.22 P2O.18P 0.065) (eg)? 2 (6.9) 
iN 

where 

N, = H,/AD,< = stability number based on the significant wave height 

Cy = empirical coefficient introduced here 

N,, = number of waves associated with the damage S starting from zero 

damage 


For depth-limited breaking waves on a sloping beach in front of the 
structure, Van der Meer (1988) proposed the use of H,., and rewrote his equation for 
relatively deep water using H,,., = 1.4H, based on the Rayleigh distribution of wave 
heights where €,, based on H, was not changed. Note that the values of H>,, were not 
tabulated in his report, so it is difficult to evaluate this shallow-water extension of his 
stability equations. In the present experiment, incident irregular wave breakers were 
spilling and plunging on the 1:20 beach slope and collapsing and surging on the 1:2 


structure slope. However, plunging waves on the beach were also observed to hit the 
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structure directly. Consequently, although it is likely that the beach slope was most 
influential on the breaking process, an appropriate slope for €,, is uncertain and Equation 
6.9 is simply adopted, although the relation does not contain €,,. The stability number 
N, in Equation 6.9 based on H, is used here because it will be easier to predict H, than 
H,,, using time-averaged surf zone models such as Battjes and Stive (1985) and Dally 
(1992) that will need to be modified to account for reflected waves (Baquerizo et al. 
1997). The empirical coefficient C, is calibrated using the 16 tests in Table 6.1. For 
eight tests with toe depth h, = 0.4 m and H, /h, = 0.23-0.39, C, = 0.97-1.21 for N,, = 
1,000 and C, = 0.99-1.19 for N,, = 3,000. For eight tests with h, = 0.2 m and H, /h, = 
0.64-0.72, Cy = 1.03-1.27 for N,, = 1,000 and C, = 1.09-1.34 for N,, = 3,000. The 
present experiment with H, /h, = 0.64-1.11 using Table 4.4 is more similar to the eight 


tests in Table 6.1 with h, = 0.2 m. So C, = 1.2 is tentatively assumed in Equation 6.9. 


Equation 6.9 is rewritten as two damage formulas in the following to 
- facilitate the representation of incident random waves in time series and spectra. The 


first damage model, based on wave time series statistics, has the form 
: b 
SS GIy | == 6.1 
| T | (6.10) 
where t = T,, N,, is the test duration for constant wave conditions. The empirical 
coefficient a, is related to cot 0, P, and C,, and b in Equation 6.10 is introduced for long 


duration tests, where b = 0.5 in Equation 6.9. The second model is based on the wave 


frequency domain statistics and has the form 
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b 
7] (6.11) 


In Equation 6.11, N,,, = H,,. /AD,,59 is the stability number based on the zeroth moment 
wave height, the test duration ¢ for constant wave conditions is normalized by the peak 
period, and a, and b are again empirical coefficients that will be obtained from data. It 
is noted that Equations 6.10 and 6.11 assume S = 0 at t= 0 and constant wave 


conditions. 


The expression for a, in Equation 6.10 can be found from Equation 6.9 
where P = 0.4 for conventional rubble-mound breakwaters. For Cy = 1.2, P=0.4, and 
cot 8 = 2, a, = 0.003 for Equation 6.10 if b = 0.5 is applicable to the present long 
duration tests. Van der Meer (1988) analyzed the five long duration tests with N,, up to 
15,000 conducted by Thompson and Shuttler (1976) and obtained the term $'/N,,°° in 
Equation 6.9 for N,, = 1,000-8,500. His data analysis indicated b < 0.5 for N,, > 8,500. 
In the following, the values of a, and b are calibrated for this experiment where N,, = 


60,000 for Series A’ and N,, = 18,000 for Series B’ and C'. 


The empirical relationship in Equation 6.10 is limited to damage due to 
incident irregular waves with constant H, and T,, starting from S=0atr=0. Asa 
result, these formulas and other existing formulas are intended for prediction of damage 


during the peak of a design storm. In order to develop an empirical procedure for more 
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realistic conditions of H, and T,, varying with time, the rate of damage increase d$'/dt is 


obtained from Equation 6.10 as 


= STarbNewieeaiva! (6.12) 
t 


which is assumed to be valid at arbitrary time t. To apply Equation 6.12 for H, and T,, 
varying with time 1, the relation must be integrated numerically with an arbitrary initial 
value of .§. For practical applications, the values of H, and T,, may be assumed to be 
constant for a short duration. Integration of Equation 6.12 for the duration of constant 


H, and T,, lasting t = t, to t=1,,, yields the mean damage at arbitrary time ¢ 
Si) = Sit) + a,NOT,'(°-1,) fort, <t<t,, (6.13) 


where S(t) = known damage at t=1,. Use of Equation 6.13 for each interval of 
constant H, and T,, in sequence allows one to compute §(1,,) for incident wave 
conditions represented by H, and T,, varying with time. Alternatively, representing the 
irregular waves at the toe of the slope using spectral parameters as shown in Equation 


6.11 yields the equation 


SO) = 0G) + OME Clan) fer a. 88St,, (6.14) 


pmo 


Equations 6.13 and 6.14 are valid for an undamaged structure (S = 0 at t= 0) or 
damaged structure, S(t)) = S; at t= 7%). Using the coefficients predicted by Van der 


Meer’s data, a, = 0.003 and b = 0.5, Equation 6.13 follows the general trend of the data 
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but underpredicts damage up to N,,=5000. For a greater number of waves, this equation 
overpredicts damage. 

The empirical coefficients in Equations 6.13 and 6.14 are recalculated in the 
following using the present experimental data, which were obtained from tests much 
longer in duration than those listed in Table 6.1. The wave heights and periods for the 
present experiment were listed in Tables 4.3 and 4.4 and the experimental data are 
summarized in Appendix B. In order to get a more gradual increase in S(t) it was 
necessary to raise a and lower b. The final coefficients were a, = 0.025 and b = 0.25 for 
Equation 6.13 and a, = 0.022 and b = 0.25 for Equation 6.14 yielding the general 
predictive relations for mean damage progression due to breaking waves for time and 


frequency wave Statistics, respectively, as 


BO 2G) > ODSQPIR CF 26") meester, (6.15) 


0.25 


SG) SSG) o OOM SI Ck a=) for L ses ty (6.16) 


n 


Figure 6.7 shows Equation 6.15 (dashed line) and Equation 6.16 (solid line) plotted 
against the damage data of Figure 5.10 for Series A’. The mean damage is well 
predicted by Equations 6.15 and 6.16 with a correlation coefficient of r = 0.99 for both 
equations for this series. Similarly, Figures 6.8 and 6.9 show that the equations predict 
damage well for Series B’ and C’, with correlation coefficients of r = 0.99 for both. 


Equations 6.15 and 6.16 are therefore generalized for Series A’, B’, and C’. The 


124 


agreement for Equations 6.15 and 6.16 is good for the three series except that the 
formulas do not predict damage stabilization adequately. This deficiency becomes 
much worse for b= 0.5. The calibrated value of b = 0.25 indicates that damage varies 
less with time ¢ in the long duration experiments using uniform-sized stone armor. The 
calibrated values of a, = 0.025 and a, = 0.022 for b = 0.25 are much larger than a, = 
0.003 for b = 0.5 estimated using data from Table 6.1. This difference is mainly caused 
by the different value of b. As an example, for t/T,, = 10,000, (¢/T,,)’ = 10 or 100 for 

b =0.25 or 0.5. Consequently, the values of a, and a, are coupled with the calibrated 


value of b. 
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Figure 6.7. Damage prediction relations compared to data for Series A’ 
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Figure 6.9. Damage prediction relations compared to data for Series Cs 


126 


Chapter 7 
EFFECT OF WAVE PERIOD AND 
ARMOR GRADATION ON DAMAGE PROGRESSION 

7.1 Experimental Setup and Test Conditions 

It is expected that the empirical coefficients in Equations 6.15 and 6.16 are, 
to some extent, a function of wave period. In addition, only uniform sized armor was 
utilized in Series A’, B’, and C’. Therefore, four test series (D’, E’, F’, and G’) were 
added to address the effect of wave period and stone gradation on damage progression. 
These additional test series até summarized in Table 7.1 along with Series A’, B’, and C’ 
for comparison. The spectral and time series incident wave parameters for all series are 
listed in Table 7.2. Damage measurements are summarized for Series D', E’, F'’, and G’ 


- in Tables B4 through B7 in Appendix B. 


Table 7.1 
Series Type Type Order (hr) 
[7 [Storm Ordering | _Uniform | High-Low | 9.0_—+| 


|G’ | __Gradation ___|__Riprap_ | tow-High | 85 
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Series D' and E' were similar to Series B’ except with different target peak 
periods. The different spectra also produced different characteristic wave heights at the 
toe of the structure. Series F’ and G’ were similar to Series B’ and D’ except that riprap 
was used as armoring. Series F’ had the same wave peak periods as Series B’ while 
Series G’ had the same wave peak periods as Series D'. Therefore, Series B’, D’, E’, F’, 
and G’ provide comparative damage measurements with systematic variation of armor 


gradation and wave period. 


As shown in Table 7.1, the order of water levels, individual storm durations, 
and series durations for Series D', E’, F’, and G’ were identical to Series B’. Four storms 
were run within each test series following 30 min of shakedown waves, with the low- 
water-level storms run first followed by the high-water-level ones. The duration of each 
of 4 storms was 2 hours yielding 8.5 hours of waves for each series. Wave heights were 
incrementally raised at each water level to simulate increasing storm severity. The 

" waves were run in 15 min bursts, with the water completely settling between bursts. As 
was done for Series A’, B’, and C’, the 2 structures were profiled every 30 min using 
eight profile rods across the center portion of each structure. Each storm series was 
repeated to yield 32 profiles alongshore for the undamaged underlayer and armor layer 
and after every 30 min of waves. As discussed in Chapter 5, the mean and standard 


deviation of E, C, L, and S for these 32 profiles were computed. 


The armor stone for Series F’ and G' was widely graded riprap with a 


median mass M,, = 256 g, nominal diameter D,,<) = (M,,/p,)'” = 4.58 cm, stone density 
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Pp, = 2.66 g/cm’, and D,,/D,; = 1.53. The mass distribution, plotted as percent finer by 
mass, is shown in Figure 7.1. The median mass was computed from a 200 stone sample 
as the value where half the total mass in the sample was greater. For comparison, the 
median mass by count for a 200 stone sample was M,,. = 128 g (i.e. half the stones 
weighed less). This latter value was not used in this thesis, although it is sometimes 
used in the field. The riprap followed the widest recommendation of the SPM (1984) of 
approximately 0.125M..-< M<4M.». For all series, the underlayer or filter layer had a 
gradation of D,./D,; = 1.32 and was sized such that (Ms) armor / (Mso) giter = 25 and 


(D 50) armor / (D 50) filter = 2. 
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Figure 7.1. Stone mass distributions for mprap and underlayer 
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7.2 Probability Density Functions of Damage, Eroded Depth, and Cover Depth 
The probability density functions of $*, E*, and C* are shown in Figures 
7.2, 7.3, and 7.4, respectively, with Series D’, E’, F’ and G’ combined onto each figure. 
These figures can be compared to Figures 5.2, 5.3, and 5.4 for Series A’, B’, and C’, 
respectively. The figures indicate that the variability of the eroded profile and eroded 
area for the new series are similar to those of Series A’, B’, and C’. The limits of 
variability remain as prescribed by Equation 5.10 or approximately in the range from 
-3.0 to 3.0. It is noted that the damage and eroded depth distributions are negatively 
skewed while the cover depth is positively skewed. This indicates that more cross 
sections have less normalized damage and larger normalized cover layer thickness. 
Relatedly, the negative tails of the damage and eroded depth distributions are less 
populated for Series D', E’, F' and G’ than for Series A’, B’, and C’, particularly below 
normalized values of -2. As will be described in the following sections, wave energy 
and therefore damage levels were lower for the new series than for Series A', B’, and C’. 
The damage increased rapidly at low damage levels, = 0-1, resulting in few damage 


measurements at these low damage levels. 
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e SeriesD' =* SeriesE' «+ SeriesF' + Series G’' —— Normal 


Figure 7.2. Probability density function for normalized damage S* for Series D’, E’, F’ 
and G' 


e SerilesD' © SeriesE' 4 SeriesF' + Series G’ — Normal 


Figure 7.3. Probability density function for normalized eroded depth E* for Series D’, 
E’, F' and G' 
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e SeriesD' > SeriesE' « SeriesF' - Series G’ —— Normal 


Figure 7.4. Probability density function for normalized cover depth C* for Series D’, E’, 
F' and G' 
7.3 Damage Variability 

Figure 7.5 shows 6, as a function of mean damage S'for Series D’, E’, F’, 
-and G' along with Equation 6.1. Series B’ is included for comparison. This figure 
indicates that Equation 6.1 underpredicts damage variability for the new series in the 
range S'= 1-4. But above and below this range the results are inconclusive. The greater 
variability in damage for the wider gradation (Series F’ and G’) is expected; but the 
reason for the greater damage variability for shorter wave periods is not clear. It appears 
there are insufficient data to make a definitive modification to Equation 6.1. Even 


though the data are slightly underpredicted by Equation 6.1, Equations 5.10 and 6.1 
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provide a description of the alongshore variability of damage that is reasonable for 


experimental results described herein. 
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Figure 7.5. Prediction of damage variability, characterized by the standard deviation, as 
a function of mean damage for Series B’, D’, E’, F’, and G’ 
7.4 Eroded Profile Prediction 

Figure 7.6 shows the mean eroded depth F for the four new series as a 
function of S along with Series B’. Equation 6.2, for prediction of £ based on Series A’, 
B', and C’, is also plotted in Figure 7.6. It is clear that this equation describes the new 
data well without modification. The standard deviation for eroded depth o; as a 
function of Sis shown in Figure 7.7 for the four new series and B’. Equation 6.3 for 6; 


as a function of S, also plotted in Figure 7.7, underestimates the variability of eroded 


134 


depth in the new data sets. This reinforces the fact that damage variability was greater 
in Series D', E’, F’, and G' than in Series A’, B’, and C’. Also, Series D’ produced more 
outlying points in this Figure than the other series. This is consistent with previous 
figures and indicates that further testing is required to better quantify the coefficient of 
variation of the damage and profile parameters, particularly for low damage levels. The 
mean normalized eroded length Z is plotted as a function of S'for data from the four 
new series in Figure 7.8. Equation 6.4 for Z as a function of Sis shown to provide an 


excellent fit to the new data. 
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Figure 7.6. Prediction of mean maximum eroded depth as a function of mean damage 
for Series B', D’, E’, F', and G' 
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Figure 7.7. Prediction of standard deviation of maximum eroded depth as a function of 
mean damage for Series B’, D', E’, F’, and G' 
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Figure 7.8. Prediction of mean maximum eroded length as a function of mean damage 
for Series B’, D’, E’, F’, and G’ 
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The mean cover depth as a function of S' was shown to be described by 
Equation 6.5. This relation, rearranged to be eroded cover depth, also provides a good 
fit to the new data, as shown in Figure 7.9. In addition, the standard deviation of 
minimum cover depth 6, was shown to be described by Equation 6.6. This equation, 
along with the new data, are shown in Figure 7.10. Equation 6.6 underpredicts the 
variability in cover depth, as expected. Again, Series D’ produced most of the outlying 


points. Data from the other series are reasonably well predicted by Equation 6.6. 


The profile and variability relationships given by Equations 5.10 and 6.1 
through 6.6, derived from Series A’, B’, and C’ data, fit the new data remarkably well, 
considering the ranges of stone grading and peak wave period in this experiment. 
Equation 6.3 and 6.6 underpredict the variability and will likely require modification as 
more data are acquired. It appears that the coefficient of variation is more scattered for 
low damage levels. The relations given above allow prediction of profile shape and 
- alongshore variability of the profile. The above relations were developed using only 


profile data without regard to incident wave conditions or water levels. 
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Figure 7.9. Prediction of mean minimum cover depth as a function of mean damage for 
Series B’, D’, E', F', and G' 
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Figure 7.10. Prediction of standard deviation of minimum cover depth as a function of 
mean damage for Series B’, D’, E’, F’, and G’ 
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7.5 Temporal Damage Development 

In Chapter 6, the generalized Equations 6.15 and 6.16 for predicting damage 
progression were highly correlated with damage data from Series A’, B’, and C’. Figures 
7.11 through 7.14 show Equations 6.15 and 6.16 plotted against data from Series D’, E’, 
F’, and G’. To recap, Series D’ and E' were nearly identical to Series B' except that the 
peak wave period was changed for the three series. Series F' and G’ were again similar 
except that the uniform armor was replaced with riprap. Two different peak wave 
periods were tested in Series F’ and G’. It can be seen that the damage progression 
equations predict the overall trend of damage progression for Series D' through G’ using 
a, = 0.025, a, = 0.022, and b = 0.25, although there are noted discrepancies. For 
example, it can be seen that damage initiation is underpredicted in all series. Equations 
6.15 and 6.16 significantly underpredict damage initiation, if only 1 or 2 stones are 
displaced at the beginning of each test series. This underprediction appears to be 
produced by the variability in damage initiation. Figures 7.15 through 7.18 show 
Equation 6.15 starting from two different points: from N,,=0 and from the first 
measured damage point past N,,= 1000. As can be seen, the prediction is much better 
when the damage is predicted after initial profile adjustment. Figures 7.11 through 7.18 
indicate that the empirical coefficients in Equations 6.15 and 6.16 may vary somewhat 
with wave period and stone gradation for the range of experimental conditions described 
herein. The initial profile adjustment may need to be accounted for in test series with 


relatively small cumulative damage. 
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Figure 7.12. Damage prediction relations compared to data for Series E’ 
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Figure 7.16. Damage prediction relations compared to data for Series E’ including 
damage initiation adjustment 
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Figure 7.17. Damage prediction relations compared to data for Series F’ including 
damage initiation adjustment 
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Chapter 8 


SUMMARY AND CONCLUSIONS 


This thesis addresses depth-limited breaking wave damage on rubble mound 
breakwaters. Chapters 1-3 are introductory. Chapter 1 introduces the topic of rubble 
mound structures and their failure modes. It is noted that deterioration models are 
required for risk and reliability analyses as part of rehabilitation studies for these 


structures. 


Chapter 2 discusses a new experiment conducted to investigate incipient 
motion of stone armor units in depth-limited breaking waves. Observations were made 
of incipient motion on both stone and sphere armor layers. Velocity measurements were 
made during these experiments both inside and just outside the armor layer using a laser 
Doppler velocimeter. It is shown that the maximum vertical convective acceleration 
across the armor layer is proportional to the square of the maximum vertical velocity. 
Equation 2.11, an incipient motion criterion similar to the Shields criterion for sediment 
motion, is derived based on Morison forcing for the dominant mode of motion: vertical 
lift under the steep breaking-wave face. Equation 2.11 is restated in terms of the 


stability number as Equation 2.12 as follows 
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<(C a, Cy 


The experimental observations and measurements are used to validate this incipient 


motion criterion for breakwater armor. 


In Chapter 3, historical damage experimental techniques, measurements, 
and models are described for rubble mound breakwaters. It is noted that there are 
inconsistencies in both historical and recent modeling techniques which have a large 
effect on the conclusions drawn. In particular, different methods of computing eroded 
area or volume produce results that differ by nearly a factor of 2. Several authors have 
made note of the fact that irregular wave experiments may not produce an equilibrium 
profile prior to failure. This conclusion was verified within this study. Many 
breakwater stability experiments have been conducted; but nearly all tests were 
conducted with nonbreaking waves, started with an undamaged structure, and 


progressed to relatively low levels of damage. 


Chapters 4-7 discuss a new experiment to investigate damage development 
on a breakwater exposed to a series of depth limited storms of varying wave energy and 
in varying water levels. It is noted that this study was motivated by the fact that no 
general irregular-breaking-wave damage progression experiments have been conducted. 
Therefore, no generally applicable relations exist that would permit prediction of 


damage progression in breaking waves. A methodology is developed for measuring 
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damage in a small-scale physical model. The methodology uses damage defined by 


Equation 3.1 as the eroded area normalized by the square of the nominal stone diameter 


In addition, three profile descriptors of engineering interest are suggested as 
follows. The minimum cover depth (Equation 5.1), maximum eroded depth (Equation 


5.2), and maximum eroded length (Equation 5.5) along the profile are normalized as 


D,50 


Data were obtained from an automated profiler on a traditional rubble mound 
breakwater section. The profiler is shown to provide improved damage measurements, 
as the cross-shore spatial sampling interval was less than 1 mm, which is considerably 
smaller than traditional sounding measurements. In addition, the profiler was automated 


making profiling very efficient. Because profiling did not require draining the wave 
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flume between tests, more tests could be accomplished and with more frequent 


profiling. 


In the experiment, damage development was measured on two identical 
breakwaters with seaward slopes of 1:2 exposed to depth-limited breaking waves on a 
1:20 beach slope. The armor was limited to a traditional two-layer thickness of 
uniformly sized stone and a widely graded riprap. The waves were run as sequences of 
storms with varying irregular wave conditions and water levels. Seven unique long- 
term damage test series were conducted. The first, Series A’, was a long duration test 
with varying water level and wave height. In this first series, the water level was 
increased midway through the test. The wave heights were increased in five increments. 
Damage was allowed to reach an observed apparent equilibrium level before the wave 
height or water level was changed. The second and third series (B’ and C’) were 
designed to yield insight into damage progression for variations in storm ordering. 
~ During Series B’, the water level was increased midway through the test, while for 
Series C’ the water level was decreased midway through. For Series B’ and C', wave 
heights were increased in four increments, each after about 4,200 waves or 
approximately 8 to 10 hr of prototype storm. Series D’ and E’ were similar to Series B’ 
except with different peak wave periods. Series F’ and G' were similar to Series B’ and 


D’' except that riprap was used as armoring. 


It was observed during the series that the damage never really stabilized or 


reached an equilibrium. Also, Series B’ and C’ indicated that the final damage was 
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similar for both series consisting of different sequences of storms of similar cumulative 
wave action. The damage data showed partial healing during damage progression, 
where eroded pockets from displaced stones were filled in by minor shifting of 


surrounding stones. 


The probability density functions of damage, eroded depth, and cover depth 
were analyzed in Chapter 5. For each series, at each point in time, the probability 


density functions of normalized values, given by Equation 5.6 as 


are shown to have ranges given by Equation 5.10 as 
DY SSS 3 | Ph S/B"S 3 a SCS 2B 


All normalized variables are in the range from -3 to 3. It is shown that series with lower 
overall damage levels had negatively skewed distributions with fewer points in the 
negative tail in the range -3 < S* < -2 due to the fact that there were few measurements 


at very low damage in the range S=0 - 1. 


Equation 6.1 was obtained through statistical analysis of the profile data to 


describe the alongshore variability of damage as 


Gy = OSS" 
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This equation, along with Equation 5.10, can be used to predict the range of damage for 
a given level of mean damage. These relations indicated that alongshore damage 
variability was greater than previously suspected and could be significant in a life-cycle 


cost analysis. 


The damaged profile shape, as characterized by the mean eroded depth F 
and mean cover depth C, was shown to be well described in terms of the mean damage 


by Equations 6.2 and 6.5 


where C, is the mean cover depth when S=0. The mean eroded length was shown to 


be described in terms of the mean damage by Equation 6.4 as 


L = 4.4/8 


The variability of eroded depth and cover depth was shown to be described by 


Equations 6.3 and 6.6, respectively, as 


6, = 0.26 - 0.00007(S - 7.8)" 


Ge = OG + 0.098 - 0.002(S - 7)? 


149 


where G,, is the value of o, when S'= 0. It is noted that Equations 6.1, 6.3, and 6.6 
underpredict variability of erosion somewhat, but more data are required to better 
quantify these relations. Equations 5.10 and 6.1 through 6.6 define the damage profile 
shape and alongshore variability as a function of mean damage without regard to wave 


conditions. 


The temporal variation of mean damage is given by Equations 6.15 and 6.16 


Sq) = S(r,) + 0.025(N,)* T° (95 - 1,°) 


a JOP G- SESE os 


SQ) 2 SG) + ORO PU eC = i) fer Ragen 
for time domain and frequency domain wave statistics, respectively. These formulas 
can be used for predicting mean damage .S, where the wave height and period vary in 
steps. The formulas fit the data well for Series A’, B’, and C’, with systematic variations 
of water depth and wave height. The equations are shown to fit the trends of damage 
progression well for Series D’, E’, F’, and G’, where the peak wave period and stone 
gradation were varied. Damage initiation, where only 1 or 2 stones moved, was 
consistently underpredicted by more than a standard deviation for Series D’ through G’. 
This appears to be due to the variability in damage initiation. It is shown that the 
prediction is significantly improved if the damage progression is begun immediately 
after the initial profile adjustment or 1000 waves. The initial profile adjustment may 


need to be accounted for in test series with relatively small cumulative damage. Also, 
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the empirical parameters in these equations may need to be varied somewhat with wave 


period and armor gradation. 


The relations given in this thesis for eroded area and damaged profile 
prediction can be used to predict the future performance of a new or damaged structure. 
They can be used to predict damage for a single design event or in a simulation to 
predict the life-cycle costs. Further, because the variability in damage is given, these 
equations can be used to predict the reliability of a given design or of a damaged 
structure. The equations described above are valid for the range of wave, water level, 
and structure conditions within this experiment. Employment of a critical stability 
number for measurable damage may be necessary to reduce the number of storms 
required in a life-cycle analysis. These relations should be conservative because the 
structures were relatively impermeable, relatively low crested, and exposed to severe 
breaking waves. But they may not be conservative for more highly grouped wave 
conditions and other structure slopes. Therefore, care should be taken in extending 
these relations for conditions not tested herein. Additional experiments are being 


conducted to evaluate the generality of these equations. 
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APPENDIX A 


SUMMARY OF WAVE MEASUREMENTS 
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Figure A.1. Wave generator command signals for six wave cases of initial experiment 
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Figure A.4. Wave 1 measured time series from wave gages 
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Figure A.6. Wave 3 measured time series from wave gages 
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Figure A.7. Wave 4 measured time series from wave gages 
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Figure A.8. Wave 5 measured time series from wave gages 
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Figure A.9. Wave 6 measured time series from wave gages 
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Figure A.10. Wave 1 incident and reflected spectra from wave gages 
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Figure A.11. Wave 2 incident and reflected spectra from wave gages 
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Figure A.12. Wave 3 incident and reflected spectra from wave gages 
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Figure A.13. Wave 4 incident and reflected spectra from wave gages 
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Figure A.14. Wave 5 incident and reflected spectra from wave gages 
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Figure A.15. Wave 6 incident and reflected spectra from wave gages 
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Figure A.16. Wave 7 incident and reflected spectra from wave gages 


167 


Offshore Array 


37.5 


3 
bi 
»~ 
(3) 
o 
Ay 
n 
~ 
=} 
o 
uo) 
3) 
a 
_ 


25.0 


(zH/zW0) (3)"s 


-12.5 


0.2 


Frequency (Hz) 


Nearshore Array 


3} 
te 
~ 
oO 
o 
Ay 
n 
~ 
q 
o 
uo} 
— 
1) 
& 
— 


37.5 


———- Reflected Spectrum 


25.0 


Ve) 
nN 
nl 


(zH/2W0) (3)"S 


-12.5 


0.2 


Frequency (Hz) 


Figure A.17. Wave 8 incident and reflected spectra from wave gages 
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Figure A.18. Wave 9 incident and reflected spectra from wave gages 
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Figure A.19. Wave 10 incident and reflected spectra from wave gages 
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Figure A.20. Wave 11 incident and reflected spectra from wave gages 
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Figure A.21. Wave 12 incident and reflected spectra from wave gages 
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Figure A.22. Wave 13 incident and reflected spectra from wave gages 
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Figure A.23. Wave 14 incident and reflected spectra from wave gages 
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Figure A.24. Wave 15 incident and reflected spectra from wave gages 
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Figure A.25. Wave 16 incident and reflected spectra from wave gages 
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Figure A.26. Wave 17 incident and reflected spectra from wave gages 
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Figure A.27. Wave 18 incident and reflected spectra from wave gages 
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Figure A.28. Wave 19 incident and reflected spectra from wave gages 
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Figure A.29. Wave 20 incident and reflected spectra from wave gages 
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Figure A.30. Wave 21 incident and reflected spectra from wave gages 
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Figure A.31. Wave 22 incident and reflected spectra from wave gages 
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Figure A.32. Wave 23 incident and reflected spectra from wave gages 
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Figure A.33. Wave 24 incident and reflected spectra from wave gages 
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Figure A.34. Wave 25 incident and reflected spectra from wave gages 
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Figure A.35. Wave 26 incident and reflected spectra from wave gages 
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